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1. Introduction 

1.1. Motivation 

This study grew out of our attempt to answer the question raised by A. Edclman 
(1997): "What is the probability p n ^ that an n x n random real Gaussian matrix 
has exactly fc real eigenvalues?" In the physics literature, an ensemble of such 
random matrices is known as GinOE - Ginibre's Orthogonal Ensemble (Ginibrc 
1965). Looking into this particular problem, we have realised that no comprehensive 
solution for the probability p n ^ can be found without undertaking an in-depth study 
(Kanzieper and Akemann 2005) of the intcgrablc structure of GinOE. The results of 
our investigation are reported in the present paper. 

1.2. Main results 

For the benefit of the readers, we collect our main results into this easy to read 
subsection with pointers to the sections containing detailed derivations of each 
statement. 

1.2.1. Real part of GinOE spectrum 

(I) Probability of Exactly k Real Eigenvalues. Let 7i be an n x n ran- 
dom real matrix whose entries are statistically independent random variables picked 
from a normal distribution N(0, 1). Then, for n — k = 21 even, the probability p n ^ of 
exactly k real eigenvalues f occurring is 

Pn,k = Pn,n-2l = ■ ■ ■ ,Pl), (1.1a) 

where p n ,n is the probability p n . n = 2~™("~ 1 )/ 4 of having all n eigenvalues real 
(Edelman 1997). The universal multivariate functions Z^e^, solely determined by the 
number I of pairs of complex conjugated eigenvalues, are so-called zonal polynomials 
(Macdonald 1998) that can be written as a sum over all partitions $ A = {l\ 1 , • • • , f g 9 ) 
of the size |A| = I 

V)(pi.-.w) = (- 1 )'«E i\h{- VJ t) ■ {llh) 

\\\=i j=i ^ v 3 j 

A few first zonal polynomials are displayed in Table 1. The arguments pj's of the 
zonal polynomials are nonuniversal § 

Pj = tr ( o,L„/2j-i) Q\ (1-lc) 

f The number of complex eigenvalues n— k = 21 is always even since the complex part of the spectrum 
consists of I pairs of complex conjugated eigenvalues. 

X The notation A = (l\ 1 , • ■ • , l g 9 ) is known as the frequency representation of the partition A of the 
size |A| = £. It implies that the part £j appears Uj times so that £ = Y2j—i &j where g is the 
number of inequivalcnt parts of the partition. In particular, the partition A = (l e ) equals (1, • • ■ , 1). 

£ times 

§ The notation tr( , a ) M denotes the trace of an (a + 1) X (a + 1) matrix M such that tr(Q a ) M = 
^2j = Mjj . Also, L^J stands for the floor function. In what follows, the ceiling function [x] will be 
used as well. 
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Table 1. Explicit examples of zonal polynomials Zne\(pi, ■ ■ ■ ,pi) as defined by 
the equation (1.16). Another way to compute Z,^-, is based on the recursion 
equation (Macdonald 1998) 

e ~\ r_iy-r-l 

^(1<)(P1,- •• ,Pl) = {(■- !) ! j Pl-r " >Pr) 

r=0 r ' 

supplemented by the formal "boundary" condition ^(i0) = 1. They are also 
tabulated in the manuscript by H. Jack (1976). 

(1 1 ) Pi 

(1 2 ) pj-P2 

(1 3 ) p?-3piP2+2p 3 

(1 4 ) p\ + 8pip 3 - 6pfp 2 + 3p| - 6p 4 

(1 5 ) p\ - 10pfp 2 + 20pfp 3 + 15pipl - 30pip 4 - 20p 2 P3 + 24p 5 



for they depend on a nonuniversal matrix q. For n = 2m even, its entries are 

r-OO 

= / t/ 2 ^" 1 e» 2 erfc(yV2) (2a + 1) 

x i^ 4 7 1 a) -\-2 y 2 ) + 2 y 2 i^ +1 (-2 y 2 )] (l.ld) 
while for n = 2m + 1 odd, 

-odd _ -even / A\m-0 TO - (-ffi)- -even /•. i \ 

2a,/3 - /3 - l- 4 J (2m)! Q '™' ^ ' 

Here, the notation erfc(t/>) stands for the complementary error function, 

2 f°° _ 2 

erfc(</>) = —= I dte * , 
while £"(</>) denote the generalised Laguerre polynomials 
L ?W ) = ^-e*|i(^c-*). 

The above result (1.1a) will be derived in Section 6. 

(II) Generating Function For Probabilities p n , k . The generating function 
G n (z) for the probabilities p nj fc is 

[n/2] 

Gn(z) = Z Wn-2e = Pn,n det [1 + Z g] [n/2J x [n/2J . (1.2) 

Equation (1.2) with the g of needed parity provides us with yet another way of 
computing the entire set of p„.fe's at once! Table 2 contains a comparison of our 
analytic predictions with numeric simulations. The result (1.2) will be proven in 
Section 6. 
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(III) Probability p n , n -2 of Exactly One Pair of Complex Conjugated 
Eigenvalues. For k = n — 2, the probability function p n ^ reduces to 

roc 

Vn,n~2 = 1Vn, n / dy y e y2 exfc(yV2) L 2 n _ 2 (- 2y 2 ) . (1.3a) 
Jo 

An alternative expression reads: 

Vn,n — 2 Pn, 



[n/2\-l 



n/2j 



(1.36) 



Here, a„ = \n/2~\ — L n /2_K and Pn{4>) stands for the Legendre polynomials 



( 1)3 J*l [(i - 6 2 y] . 



The leading large-n behaviour of the probability p n ,n-2 is given by 



3«+l/2 



i IK n 



■Pn 



(1.3c) 



The above three results will be derived in Section 4.1, Section 7.1 and Section 7.2, 
respectively. 

Complex part of GinOE spectrum 

(IV) Joint Probability Density Function of All Complex Eigenvalues 
Given There Are k Real Eigenvalues. Let Hk be an n x n random real matrix 
with k real eigenvalues such that its entries are statistically independent random vari- 
ables picked from a normal distribution IM(0, 1). Then, the joint probability density 
function (j.p.d.f.) of its 2i = n — k complex eigenvalues is 



Pn k {zi,---,Zi) 



Pn, 



X 



Y\_ erfc 



i=i 



iV2 



pf 



2£x2£ 



(1.4) 



Here, pf denotes the Pfaffian. The above j.p.d.f. is supported for (Rezi, • • • , Rez^) G 
R £ , and (Imzi, • • • ,Imzt) <G (M + ) £ . The antisymmetric kernel T> n (z,z') is given 
explicitly by (3.12) - (3.18) of Section 3 where the statement (1.4) is proven. 



(V) Correlation Functions of Complex Eigenvalues in The Spectra Free 
OF Real Eigenvalues. Let Tig be an n x n random real matrix with no real eigen- 
values such that its entries are statistically independent random variables picked from 
a normal distribution N(0, 1). Then, the p-point correlation function (1 < p < £) of 
its complex eigenvalues, defined by (2.7), equals 



n 



3=0 ' 3 



X pf 



UU r(i/2) ii 

Kt(Zi,Zj) Ki(Zi,Zj) 
Ke(z t ,Zj) K e (z t ,Zj) 



iV2 



exp 



(1.5a) 



2px2p 
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Table 2. Exact solution for P12 k (first and second column) compared to numeric 
simulations (third column) performed by direct diagonalisation of 1, 000, 000 of 
12 X 12 matrices. 



Analytic solution 



Exact 



Approximate 



Numeric 
simulation 



10 



12 



29930323227453-20772686238032^2 
17592186044416 

3(1899624551312^2-2060941421503) 
4398046511104 

3(2079282320189-505722262348^2) 
8796093022208 

252911550974^-27511352125 
4398046511104 

15(1834091507- 10083960^2) 
17592186044416 

3(1260495^2-512) 
2199023255552 



8589934592 



0.031452 
0.426689 
0.465235 
0.075070 
0.001552 

0.000002 
0.000000 



0.031683 
0.427670 
0.464098 
0.075021 
0.001526 
0.000002 
0.000000 



Here, n = 2£ and the 'pre-kernel' Kg equals 

e-i 

P2j (Z)p2j+1 (V) - p 2 j (z')p 2 j+l( z ) 



e-i 1 

Ke(z,z') = 



3=0 



(1.5b) 



The polynomials Pj(z) are skew orthogonal in the complex half-plane Imz > 0, 

(P2j+i,P2k)c = -(P2k,P2j+i)c =ir 3 S :jk , (1.5c) 

(P2j+l,P2k+l)c = (P2j,P2k)c = 0, (1.5(f) 

with respect to the skew product 

' ' z — z" 



iV2 



exp 



z 2 + z 2 



[f(z)g(z)-f(z)9(z)). (1.5e) 



(f,9)c = / d 2 zerk 

For detailed derivation, a reader is referred to Section 8 which also addresses the prob- 
lem of calculating the probability p„ ;U to find no real eigenvalues in the spectrum of 
GinOE. 



1.2.3. How to integrate a Pfaffian? 

All the results announced so far would have not been derived without a Pfaffian inte- 
gration theorem that we consider to be a major technical achievement of our study. 
Conceptually, it is based on a new, topological, interpretation of the ordered Pfaffian 
expansion as introduced in Section 5. 
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(VI) The Pfaffian Integration Theorem. Let dn(z) be any benign measure 
on z e C, and the function Q n {x, y) be an antisymmetric function of the form 

_^ n— 1 

Q n {x,y) = - ^2 <i]( x ) i^jkQkiy) (i-6o) 



j,k=0 



where q/s are arbitrary polynomials of j-th order, and (1 is an antisymmetric matrix. 
Then the integration formula 



f) n 



dn(zj) pf 



Qn{Zi, Zj) 
Qn{Zi, Zj) 



i(Zi, Zj) 
l(Zi, Zj) 



2£x2£ 



1 -1 1 -I 

2 tr (o,n-i)« , • • • , -tr( 0i „_i)t; 



(1.66) 



holds, provided the integrals in its l.h.s. exist. Here, Zni\ are zonal polynomials whose 
I arguments are determined by a matrix v with the entries 

n— 1 „ 

v a .p = i ^2 ^ a - k \ d ^( z )[lk(z)qp(z) - qp(z)q k (z)}. (1.6c) 
k=0 Jzec 

This theorem that can be viewed as a generalisation of the Dyson integration theorem 
(Dyson 1970, Mahoux and Mehta 1991) will be proven in Section 5. 

Interestingly, the Pfaffian integration theorem is not listed in the classic reference 
book (Mehta 2004) on the Random Matrix Theory (RMT). Also, we are not aware 
of any other RMT literature reporting this result which may have implications far 
beyond the scope of the present paper. 



1.3. A guide through the paper 

Having announced the main results of our study, we defer plunging into formal 
mathematical proofs of the above six statements until Section 3. Instead, in Section 
2, we deliberately draw the reader's attention to a comparative analysis of GinOE 
and two other representatives of non-Hermitean random matrix models known as 
Ginibre's Unitary (GinUE) and Ginibre's Symplectic (GinSE) Ensemble. Starting 
with the definitions of the three ensembles, we briefly discuss their diverse physical 
applications, pinpoint qualitative differences between their spectra, and present a 
detailed comparative analysis of major structural results obtained for GinUE, GinSE 
and GinOE since 1965. We took great pains to write a review-style Section 2 in order 
(i) to help the reader better appreciate a profound difference between GinOE and the 
two other non-Hermitean random matrix ensembles on both qualitative and structural 
levels as well as (ii) place our own work in a more general context. 

A formal analysis starts with Section 3 devoted to a general consideration of 
statistics of real eigenvalues. Its first part, Section 3.1, summarises previously known 
analytic results (Edelman 1997) for the probability function p n _ k of the fluctuating 
number of real eigenvalues in the spectrum of GinOE. Section 3.2 deals with the 
joint probability density function of complex eigenvalues of GinOE random matrices 
that have a given number of real eigenvalues. The Pfaffian representation (1.4) is 
the main outcome of Section 3.2. This result is further utilised in Section 3.3 where 
the probability function p n ^ is put into the form of a 'Pfaffian integral' (3.19). The 
analysis of the latter expression culminates in concluding that the Dyson integration 
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theorem, a standard tool of Random Matrix Theory, is inapplicable for treating the 
Pfafhan integral obtained. The latter task will be accomplished in Section 5. 

Section 4 attacks the probability function p n ^ for a few particular values of k. 
The probabilities p n , n -2, Pn,n-4 and p n , n -6 OI one > two, and three pairs of complex 
conjugated eigenvalues occurring are treated in Section 4.1, Section 4.2 and 4.3, 
respectively. This is done by explicit calculation of the Pfafhan in (3.19) followed 
by a term-by-term integration of the resulting Pfafhan expansion. In Section 4.4, 
we briefly discuss a faster-than-exponential growth of the number of terms in this 
expansion caused by further decrease of k. 

Section 5, devoted to the Pfaffian integration theorem, is central to the paper. 
Its main objective is to introduce a topological interpretation of the terms arising in 
a permutational expansion of the Pfaffian in the l.h.s. of (1.66). Such a topological 
interpretation turns out to be the proper language in the subsequent proof of the 
Pfaffian integration theorem. In Section 5.1, the Pfaffian integration theorem is 
formulated and discussed in the light of the Dyson integration theorem. In Section 5.2, 
an ordered permutational Pfaffian expansion is defined and interpreted in topological 
terms. The notions of strings, substrings, loop-like strings and loop-like substrings for 
certain subsets of terms arising in an ordered Pfaffian expansion are introduced and 
illustrated on simple examples in Sections 5.2.1 and 5.2.2. Further, equivalent strings 
and equivalent classes of strings are defined and counted. The issue of decomposition 
of strings into a set of loop- like substrings is also considered in detail (Lemma 5.2). 
Section 5.2.3 is devoted to the counting of loop-like strings. In Section 5.2.4, the 
notion of adjacent strings is introduced and illustrated. Adjacent strings are counted 
in Lemma 5.4. Their relation to loop-like strings is discussed in Lemma 5.5. Section 
5.2.5 deals with a characterisation of adjacent strings by their handedness; adjacent 
strings of a given handedness are counted there, too. In Section 5.2.6, equivalent 
classes of adjacent strings are defined, counted, and explicitly built. Section 5.3 and 
Section 5.4 are preparatory for Section 5.5, where the Pfaffian integration theorem is 
eventually proven. For the readers' benefit, a vocabulary of the topological terms we 
use is summarised in Table 3. 

Section 6 utilises the Pfaffian integration theorem to obtain a general solution for 
the sought probability function p n (Section 6.1), derive a detcrminantal expression 
for the entire generating function of p n ,kS (Section 6.2), and address the issue of integer 
moments of the fluctuating number of real eigenvalues in GinOE spectra (Section 6.3). 

Section 7 is devoted to the probability p n .n-2 for two complex conjugated pairs 
of eigenvalues to occur; a large-n analysis of p n , n -2 is also presented there. 

Section 8 discusses the Pfaffian structure of the p-point correlation functions of 
complex eigenvalues belonging to spectra of a subclass of GinOE matrices without real 
eigenvalues. The same section addresses the problem of calculating the probability p n $ 
to find no real eigenvalues in spectra of GinOE. 

Section 9 contains conclusions with the emphasis placed on open problems. The 
most involved technical calculations are collected in four appendices. 

2. Comparative analysis of GinOE, GinUE, and GinSE 

2.1. Definition and consequences of violated Hermiticity 

Ginibre's three random matrix models - GinOE, GinUE, and GinSE - were derived 
from the celebrated Gaussian Orthogonal (GOE), Gaussian Unitary (GUE), and 
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Gaussian Symplectic (GSE) random matrix ensembles in a purely formal way by 
dropping the Hcrmiticity constraint. Consequently, the non-Hcrmitean descendants 
of GOE, GUE and GSE share the same Gaussian probability density function 



for a matrix "Hamiltonian" H. £ T^n) to occur; the constant ap is chosen to be 
0/3 = 2 — (5/3,2 (this slightly differs from the convention used in the original paper by 
Ginibrc). However, the spaces Yp on which the matrices H. vary are different: Ti(n), 
T2(n), and Y^iri) span all n x n matrices with real (GinOE, j3 = 1), complex (GinUE, 
(3 = 2), and real quaternion (GinSE, (3 = 4) entries, respectively. 

The violated Hermiticity, 1~C ^ 7"£, brings about the two major phenomena: 
(i) complex-valuedness of the random matrix spectrum and (ii) splitting the random 
matrix eigenvectors into a bi-orthogonal set of left and right eigenvectors. Statistics of 
complex eigenvalues, including their joint probability density function, and statistics 
of left and right eigenvectors of a random matrix 7i e T@ drawn from (2.1) are of 
primary interest. 

Only the spectral statistics will be addressed in the present paper. (For studies 
of the eigenvector statistics in GinUE, the reader is referred to the papers by Chalkcr 
and Mchlig (1998), Mehlig and Chalker (2000), and Janik et al (1999). We are not 
aware of any results for eigenvector statistics in GinSE and GinOE.) 

2.2. Physical applications 

While the ensemble definition (2.1) was born out of pure mathematical curiosity || , non- 
Hermitean random matrices have surfaced in various fields of knowledge by E. Wigner's 
"miracle of the appropriateness" (Wigner 1960). From the physical point of view, 
non-Hcrmitean random matrices have proven to be as important as their Hermitean 
counterparts. (For a detailed exposition of physical applications of Hermitean RMT 
we refer to the review by Guhr et al (1998)). 

Random matrices drawn from GinUE appear in the description of dissipative 
quantum maps (Grobe et al 1988, Grobe and Haake 1989) and in the characterisation 
of two-dimensional random space-filling cellular structures (Le Caer and Ho 1990, Le 
Caer and Delannay 1993). 

Ginibrc's Orthogonal Ensemble of random matrices arises in the studies of 
dynamics (Sommers et al 1988, Sompolinsky et al 1988) and of synchronisation effect 
(Timme et al 2002, Timme et al 2004) in random networks; GinOE is also helpful in 
the statistical analysis of cross-hemisphere correlation matrix of the cortical electric 
activity (Kwapieh et al 2000) as well as in the understanding of inter-market financial 
correlations (Kwapieh et al 2006). 

All three Ginibre ensembles (GinOE, GinUE, GinSE) arise in the context of 
directed "quantum chaos" (Efetov 1997a, Efetov 1997b, Kolesnikov and Efetov 1999, 
Fyodorov et al 1997, Markum et al 1999). Their chiral counterparts (Stephanov 1996, 
Halasz et al 1997, Osborn 2004, Akemann 2005) help elucidate universal aspects of the 
phenomenon of spontaneous chiral symmetry breaking in quantum chromodynamics 
(QCD) with chemical potential: the presence or absence of real eigenvalues in the 
complex spectrum singles out different chiral symmetry breaking patterns. For a 
review of QCD applications of non-Hermitean random matrices with built-in chirality, 
the reader is referred to Akemann (2007). 

|| J. Ginibre, private communication. 




(2.1) 
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Figure 1. Numerically simulated distributions of complex eigenvalues in GinUE 
(left panel), GinSE (middle panel), and GinOE (right panel). Three different 
eigenvalue patterns are clearly observed. Eigenvalues are scattered almost 
uniformly in GinUE, depleted from the real axis in GinSE, and accumulated 
along the real axis in GinOE. 

Other recent findings (Zabrodin 2003) associate statistical models of non- 
Hermitcan normal random matrices with integrable structures of conformal maps 
and interface dynamics at both classical (Mineev-Weinstein et al 2000) and quantum 
(Agam et al 2002) scales. For a comprehensive review of these and other physical 
applications, the reader is referred to the survey paper by Fyodorov and Sommers 
(2003). 

2. 3. Spectral statistical properties of Ginibre 's random matrices 

A profound difference between spectral patterns of the three non-Hermitean random 
matrix models has been realised long ago. Anticipated in the early papers by Ginibre 
(1965) and Mehta and Srivastava (1966), it was further confirmed analytically by 
using varied techniques % (Edelman 1997, Efetov 1997a, Efetov 1997b, Kolesnikov 
and Efetov 1999, Kanzieper 2002a, Kanzieper 2002b, Nishigaki and Kamenev 2002, 
Splittorff and Verbaarschot 2004, Kanzieper 2005, Akemann and Basile 2007). 

Qualitatively, there is a general consensus that (i) the spectrum of GinUE is 
approximately characterised by a uniform density of complex eigenvalues. This is 
not the case for the two other ensembles, (ii) In GinSE, the density of complex 
eigenvalues is smooth but the probability density of real eigenvalues tends to zero. This 
corresponds to a depletion of the eigenvalues along the real axis, (iii) On the contrary, 
the density of eigenvalues in GinOE exhibits an accumulation of the eigenvalues along 
the real axis. It is the latter phenomenon that will be quantified in our paper. 

Our immediate goal here is to highlight the inter-relation between these quali- 
tative features of the complex spectra and the formal structures underlying Ginibre's 
random matrix ensembles. To this end, we present a brief comparative review of the 
major structural results obtained for all three Ginibre's ensembles (GinUE, GinSE, 
and GinOE) since 1965, in the order of increasing difficulty of their treatment. 

Joint Probability Density Function of All n Eigenvalues. In this sub- 
section, we collect explicit results for the joint probability density functions of all n 

K The difference in spectral patterns of non-Hermitean chiral random matrix models arising in the 
QCD context was first studied numerically by Halasz et al (1997). A review of recent theoretical 
developments can be found in Akemann (2007). 
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complex eigenvalues of a random matrix Ti e Tp(n) drawn from any of the three 
Ginibrc random matrix ensembles. 

• GinUE: The spectrum of a random matrix Ti € T2(n) drawn from GinUE 
consists of n complex eigenvalues (z\,---z n ) whose joint probability density 
function mirrors (Ginibre 1965) that of GUE (sec, e.g., Mehta 2004), 

(n \ ~~ n n 

e=i / £i>£ 2 =i e=i 
Similarly to the electrostatic model introduced by E. Wigner (1957), J. Gini- 
bre pointed out that the j.p.d.f. (2.2) can be thought of as that describing the 
distribution of the positions of charges of a two-dimensional Coulomb gas in an 
harmonic oscillator potential U{z) = \z\ 2 /2, at the inverse temperature 1/T = 2. 

• GinSE: The spectrum of a random matrix Ti € T^n) drawn from GinSE 
consists of n pairs of complex conjugated eigenvalues (z\, z\, ■ ■ ■ , z n , z n ). The 
corresponding joint probability density function was derived by Ginibre (1965), 

(n \ — n 

(27rrnin(2*-l)! U \z h - z e2 \ 2 \z £l - z e ,\ 2 

e=i / ii>i 2 =i 

n 

x JJ \zi - z e \ 2 exp(-zeze). (2.3) 
i=i 

Notice that the factor n"=i \ z e ~ m (2-3) is directly responsible for the de- 
pletion of the eigenvalues along the real axis. For GinSE, a physical analogy with 
a two-dimensional Coulomb gas is much less transparent; it has been discussed 
by P. Forrester (2005). 

• GinOE: Contrary to GinUE and GinSE, the complex spectrum (w\, • • • , w n ) of a 
random matrix W € Ti (n) drawn from GinOE generically contains a finite frac- 
tion of real eigenvalues; the remaining complex eigenvalues always form complex 
conjugated pairs. Indeed, no other option is allowed by the real secular equation 
det (w — Ti.) — determining the eigenvalues of Ti. 

This very peculiar feature of GinOE, that we call accumulation of the eigenvalues 
along the real axis, can conveniently be accommodated by dividing the entire 
space Ti(n) spanned by all real n x n matrices Ti £ Ti(n) into (n + 1) mutually 
exclusive sectors Ti(n/fc) associated with the matrices Tik C Ti having exactly k 
real eigenvalues, such that Ti(n) ~ Ufe=o Ti(n/fc). The sectors Ti(n/fc), charac- 
terised by the partial j.p.d.f. 's -pHeTi(n/fc) {wi, ■ ■ ■ , w n ), can be explored separately 
because they contribute additively to the j.p.d.f. of all n eigenvalues of Ti from 
Ti(n): 

n 

pW(wi, ■ •• ,W n ) = ^2P-HeT(n/k)(wi, ■ ■ -,w n ). (2.4) 

fe=0 

In entire generality, the partial j.p.d.f. 's have been determined by Lehmann and 
Sommers (1991) who proved, a quarter of a century after Ginibre's work, that 
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the fc-th partial j.p.d.f. (0 < k < n) equals 

^7teT(n/fe)(^l) 1 " " ) A*,; z \i " " " > z e) 

2^-n(n+l)/4 k k 

= vm n" 1 ru/2) n i A * - A ^'i II c M-^m 



x n ii( a j - ^)( a j - *o n \^ - z j\ 2 \ z i - * 



2 



x n (*j - crfc (^^J exp . (2.5) 

Here, the parameterisation (w\, ■ ■ ■ , io„) = (Ai, • • • , A^; z\, zi, ■ ■ ■ , ze, ze) was used 
to indicate that the spectrum is composed of k real and 2i complex eigenval- 
ues so that k + 21 = n. The above j.p.d.f. is supported for (Ai, • • • , A fc ) € R k , 
(Rezi, • • • ,Reze) € M £ , and (Imzi,---,Im^) e (M+) £ . Notice that since the 
complex eigenvalues come in i conjugated pairs, the identity n = 2£ + k implies 
that half of the sets T\(n/k) are empty: this happens whenever n and k are of 
different parity. 



In writing (2.5), we have used a representation due to Edelman (1997) who re- 
discovered the result by Lehmann and Sommers (1991) a few years later. The 
particular case k — n oi (2.5), corresponding to the matrices H € T(n/n) with all 
eigenvalues real, was first derived by Ginibre (1965). No physical interpretation 
of the distribution (2.5) in terms of a two-dimensional Coulomb gas is known as 
yet. 



Eigenvalue Correlation Functions and Inapplicability of the Dyson 
Integration Theorem to the Description of GinOE. Spectral statistical 
properties of random matrices can be retrieved from a set of spectral correlation 
functions defined as 

R(f)( Zl ,..., Zp .,n) = 7 -^— J] d 2 Zj pW( Zl ,---,z n ), (3 = 2,4 (2.6) 

for GinUE (/? = 2) and GinSE (/? = 4), and 
R p ,q k (Ai, • • • , A p ; zi, • • • , z q ; n) — ^ - — ^ _ ^ 

x n / dx i n / rf2 

z m -pHeT(n/fe)(Ai, ■ ■ ■ , Afe ; z\ , ■ ■ ■ , ze), (2.7) 

j=p+l •' R m=q+l , ' C 

for GinOE (/3 = 1). The GinOE correlation function refers to the spectrum of matrices 
Tik C H. having exactly k real eigenvalues. 

The analytic calculation of the above correlation functions is one of the major 
operational tasks of the non-Hermitean RMT. Whenever feasible, such a calculation 
either explicitly rests on or can eventually be traced back to the three concepts: (i) a 
determinant (or Pfaffian) representation of the j.p.d.f.'s of all eigenvalues, (ii) a projec- 
tion property of the kernel function associated with the aforementioned determinant 
representation, and (iii) the Dyson integration theorem (Dyson 1970, Mahoux and 
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Mchta 1991) that makes use of both (i) and (ii). Highly successful in the Hcrmitcan 
RMT, the above three concepts are not always at work in the non-Hermitean RMT. 
Particularly, the Dyson integration theorem, being effective for GinUE and GinSE 
(see below), fails to work for GinOE. 

It will be argued that a mixed character of the GinOE spectrum consisting of both 
complex and purely real eigenvalues is the direct cause of the failure. For the readers 
convenience as well as for the future reference, we cite the Dyson integration theorem 
below + (Mehta (1976); see also Theorem 5.1.4 in Mehta's book (2004)). 



Theorem 2.1 (Dyson integration theorem). Let f(x,y) be a function with real, 
complex or quaternion values, such that 

f(x,y) = f(y,x), (2.8a) 

where f = f if f is real, f is the complex conjugate if it is complex, and f is the dual 
°f f if H "is quaternion. Assume that 

dir(y)f(x,y)f(y,z) = f(x, z) + A f(x, z) - f(x, z) A, (2.86) 

where X is a constant quaternion and dir is a suitable measure. Let [f(xi, Xj)] n xn 
denote the nx n matrix with its (i,j) element equal to f(xi,Xj). Then, 

J dn(x n ) det[f(xi, Xj)] nxn =(c-n + 1) dct[/(x i; Zj)](„-i)x(n-i) ( 2 -8c) 
with 

c= J dn(x) f(x,x). (2.8(f) 

When f(x,y) is real or complex, the quaternion constant A vanishes. For f(x,y) 
taking quaternion values, det should be replaced by qdet ; the quaternion determinant 
(Dyson 1972). 



This theorem prompts the following definition. 



Definition 2.1. A function f(x,y) satisfying the first and the second equation in the 
Dyson integration theorem is said to obey the projection property. 



Being equipped with the above reminder, we are ready to present, and discuss, a 
collection of formulae available for the p-point correlation functions in Ginibre's en- 
sembles. 



• GinUE: The joint probability density function (2.2) of all n eigenvalues is 
reducible to a determinant form (Ginibre 1965) 



pV>{z u ■ ■ ■ , z n ) = - det K^(z k , z e ) ]J e 



(2.9) 



+ While the formulation in Mehta (2004) refers to the flat measure dn(x) = dx, the Theorem 2.1 
stays valid for any benign measure dn(x). 
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(2) 

with K n '(z, z) being a two-point scalar kernel 

Since it obeys the projection property, the Dyson Integration Theorem brings a 
determinant expression for the p-point correlation function: 



(2.10) 



( Zl ,---,z p ;n) = det \k ^ (z k , z t j\ f[ e 



pxp 



(2-11) 



These results, first derived by J. Ginibre (1965), provide a comprehensive descrip- 
tion of spectral fluctuations in GinUE. 

GinSE: The joint probability density function (2.3) of all n eigenvalues is 
reducible to a quaternion determinant form (Mehta and Srivastava 1966) 



PM( Z1 , ■■■ ,z n ) = — qdet \K^(z k , z e )} jj^ " ^ e ~ 

n\ L J nxn -"- - 1 - 



(2.12) 



where K n ^ {z, z') is a quaternion whose 2x2 matrix representation reads 
(Kanzicpcr 2002a) 

K^(Z,Z>) K%\Z,Z') r 



e[K^(z,z')} 



(2.13) 



Here, 



n— 1 k 



y 2k+l(„t\2t 



(z') 2£ - {z') 



(2fc + l)!!(2f)!! 



(2.14) 



k=0 £=0 

Alternatively, but equivalently, (2.12) can be reduced to the Pfaffian form 
(Akemann and Basile 2007) 

(4)/ \ (4)/ — \ 

Kn \Zii Zj) K n \Z%^ Zj) 

Kn (Zii Zj) ^n {Zi : Zj) 



1 



PW(z 1; ...,z„) =-pf 



J 2nx2n 



which is instructive to compare with (1.4). 



(2.15) 



As soon as the quaternion kernel K n 4 ' (z, z') satisfies the projection property, 
the p-point correlation functions take a quaternion determinant /Pfaffian form: 

p 

R p A \z u ■■■,z p ;n) = qdet K^(z k , z t ) Y[(zj - Zj) e~ z ^ . (2.16) 

L J pxp 

This result is due to Mehta and Srivastava (1966). 



GinOE: To the best of our knowledge, structural aspects of correlation functions 
in GinOE have never been addressed (see, however, a recent paper by Sinclair 
(2006)); consequently, no analogues of the above GinUE and GinSE formulae 
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[(2.11) and (2.16)] are available. This gap will partially be filled in the present 
paper where we derive a quaternion determinant /Pfaffian expression (1.4) for 
the j.p.d.f. of all complex eigenvalues of a random matrix Ti e T\{n/k). 
Importantly, the 2x2 kernel therein does not possess the projection property, 
hereby making the Dyson integration theorem inapplicable for the calculation 
of associated correlation functions. We reiterate that a mixed character of the 
GinOE spectrum, composed of both complex and purely real eigenvalues is behind 
the statement made *. 



Quantifying the Qualitative Differences Between Spectra of GinUE, 
GinSE, and GinOE. The mean density of eigenvalues R[^\z = x + iy;n) is the sim- 
plest spectral statistics exemplifying differences in spectral patterns of GinUE, GinSE, 
and 

GinOE. Below we collect, and comment on, the exact and the largc-n results for 
the mean density of eigenvalues of Ti. £ Tg(n). 



• GinUE: In accordance with (2.11) and (2.10), the exact result for the mean 
spectral density reads (Ginibre 1965) 

^.ro^+jfl (21?) 

Here, r(a, <j>) is the upper incomplete gamma function 

T(a,(/>) = / dtt^e" 1 . 
J<p 

For a large- n GinUE matrix and z = x + iy fixed jj, the mean spectral density 
approaches the constant 

4 2) (z;n> 1) ~ -. (2.18) 

This suggests that that the eigenvalues of a large-n GinUE matrix are distributed 
almost uniformly within the two-dimensional disk of the radius ^n. More 
rigorously, this statement follows from the macroscopic limit of (2.17) 

Y™R?\z = z^n)= 1 -{ I *W<\ (2.19) 

known as the Girko circular law (Girko 1984, Girko 1986, Bai 1997). The density 
of eigenvalues away from the disk is exponentially suppressed (Kanzieper 2005). 



• GinSE: In accordance with (2.13) - (2.16), the exact result for the mean spectral 
density reads (Mehta and Srivastava 1966, Kanzieper 2002a) 

4 4) (z; n) = 2y e~ (x2+y2) Im ^ (x + iy,x- iy) (2.20) 

* Intriguingly, it will be shown in Section 8 that the matrices H S Ti(n/0) exhibit GinSE-likc 

correlations; this contrasts the well known correlations of the GOE type (Ginibre 1965) for the 
matrices H S Ti (n/n). 
(j That is, z does not scale with n 3> 1. 
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A H ^ 



Figure 2. Profiles of eigenlevel densities (z;n), plotted as functions of the 
complex energy z = x+iy for n fixed, show a nearly uniform eigenlevel distribution 
in GinUE (left panel), a depletion of eigenvalues along the real axis in GinSE as 
exemplified by the density drop at y = (middle panel), and accumulation of real 
eigenvalues in GinOE displayed as a wall at y = (right panel); the wall height 
imitates the density of real eigenvalues. 

with given by (2.14). For a large-n GinSE matrix and z = x + iy fixed, it 
reduces to (Kanzieper, 2002a) 

4 4) {z; n » 1) ~ -j= e- 2y2 erfi(y V2) . (2.21) 
Here, erfi (</>) is the imaginary error function 

2 2 

erfi(0) = -= / die 1 . 
V 71 " Jo 

Both results [(2.20) and (2.21)] suggest that the mean eigenvalue density is no 
longer uniform but exhibits a depletion of eigenvalues along the real axis. Similarly 
to the GinUE, the mean spectral density in GinSE is suppressed away from a disk 
of the radius ^/2n as suggested by the circular law 

hm R[ 4 \zV2n-;n) = ±{ J IJ |!j < ] (2.22) 
«^oo 1 v ' ; 2tt \ if \z\ > 1 v ; 

due to Khoruzhenko and Mezzadri (2005). 

• GinOE: A mixed character of the GinOE spectrum consisting of both complex 
and real eigenvalues makes the RMT techniques based on the Dyson integration 
theorem inapplicable to the description of spectral statistical properties of GinOE. 
To evaluate the mean eigenvalue density in the finite-n GinOE, a totally different 
approach has been invented by A. Edclman and co-workers (Edelman et al 1994, 
Edclman 1997). Starting directly with the definition (2.1) taken at /? = 1 and 
applying the methods of multivariate statistical analysis (Muirhead 1982), these 
authors have separately determined the exact mean densities f of purely real 
eigenvalues (Edelman et al 1994) 

f The two are related to the (p, q) correlation functions Rp, q k for the matrices TH.^ restricted to have 
exactly k real eigenvalues [see the definition (2.7)] as follows: 

n L«/2J 

<Li(^) = %)E R M fc) (^)' M!L mP iex(^«)= E <? n - 2f) (*;n)- 

fc=l i=\ 
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l.realV 



r(n-l,.x 2 ) 2»/ 2 - 3 / 2 , ,„_! _ j2/2 (n-l x 



2 * 



r(n- 1) r(n- 1) 1 1 ' v 2 ' 2 y 

(2.23) 

and of strictly nonreal J eigenvalues (Edelman 1997): 

< omplcx (^) = f- ^^^ Ve^eMv^)- (2-24) 

Here x and y are real and imaginary parts of z = x + iy. The function 7(0, 4>) in 
(2.23) is the lower incomplete gamma function 



7(a,0) = / dtt a ~ L e 
Jo 

Importantly, no reference was made to the j.p.d.f. (2.4) and (2.5) in deriving 
(2.23) and (2.24). 

For a large-n GinOE matrix and z = x + iy fixed, the above two formulae yield 
the mean density of eigenvalues in the form 

R { l\z-n » 1) ~ -L<J(I/) + ye^ 2 crfc (yy/2). (2.25) 

V27T V 7T 

Similarly to GinUE and GinSE, the circular law (Girko 1984, Sommers et al 1988, 
Bai 1997, Edelman 1997) 

Jim R^V-n;n)= l -{ \ * |!| < \ (2.26) 

holds. 

2.4- Statistical description of the eigenvalue accumulation in GinOE 

In essence, the approach culminating in the explicit formula (2.23) for the mean density 
of real eigenvalues represents the simplest possible quantitative description of the 
phenomenon of eigenvalue accumulation along the real axis. Complcmcntarily, the 
one might be interested in the full statistics of the number Af r (Af c ) of real (complex) 
eigenvalues occurring in the GinOE spectrum. In the latter context, the result (2.23) 
can only supply the first moment of Af r - the expected number E n — E[A/" r ] of real 
eigenvalues. Indeed, integrating out the eigenlevel densities (2.23) and/or (2.24) over 
the entire complex plane, 

E n = d 2 zR ( ll c Jz;n) = n- / d 2 z ^(z; n), 

JC JC\R 

Edelman et al (1994) have obtained the remarkable result 

E n ^l + V2 2F ^\ 1 / 2 ;^ (2.27) 
2 B(n,l/2) 1 ' 

expressed in terms of the Gauss hypergeometric function and the Euler Beta function. 
As n — > 00, it furnishes the asymptotic series 

[2n / 3 3 27 499 ml *\ 1 

^ - V V 1 - ^ - 1^ + 1024^ + 32W + ° {n } + 2- (2 - 28) 



a-l„-t 



X In the formulae, we use the subscript "complex" to identify eigenvalues with zero real part. 
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The nonperturbative -y/^-clependence of leading term in (2.28) had been earlier 
detected numerically by Sommers et al (1988). 

What is about a more detailed statistical description of the number N r (-A/" c ) 01 
real (complex) eigenvalues? Can all the moments and the entire probability function 
of the discrete random variable Af r (Af c ) be determined? The preceding discussion, 
particularly highlighting a patchy knowledge of the spectral correlations in GinOE, 
suggests that both spectral characteristics, which can be expressed as § 

W] = V S(p, k) f d 2 Zl ■■■ ( d 2 z k 4 1} rcal (zi, ■ • • , z k ; n) (2.29) 

and 

E[A/?]=£;S(p,*) / d 2 Zl ---( d 2 z e R { ll ompleK (z u ---,z f ,n), (2.30) 

are not within immediate reach. First, only one correlation function out of p involved 
in either (2.29) or (2.30) is explicitly known. Second, even if all the multipoint 
correlation functions were readily available, the integrations in (2.29) and (2.30) would 
not be trivial either because of the anticipated inapplicability of the Dyson integration 
theorem as discussed in Section 2.3. 

Some of the difficulties outlined here will be overcome in the present paper. For 
a list of our major results the reader is referred back to the Section 1.2. Proofs and 
derivations are given in the Sections to follow. 



3. Statistics of real eigenvalues in GinOE spectra: A PfafRan integral 
representation for the probability p n ^ 

3.1. Generalities and known results 

Instead of targeting the moments E[A/^] of a fluctuating number of real eigenvalues 
as discussed in the previous section, we are going to directly determine the entire 
probability function p n _ k = Prob(7V r = k). The definition of the p n ,k, describing 
the probability of finding exactly k real eigenvalues in the spectrum of an n x n real 
Gaussian random matrix, can be deduced from (2.4) and (2.5), 

p nM = Prob(?i G T(n/k)) = JJ / d\ t JJ / d 2 z 

jPneT( n /k)- (3-1) 

Here, £ is the number of pairs of complex conjugated eigenvalues in the spectrum of 
an n x n matrix 7i g T(n/k) having exactly k real eigenvalues, and P-HeT(n/k) is the 
j.p.d.f. of all n eigenvalues of such a matrix. Obviously, the identity n — k + 2£ holds. 

Previous attempts to determine the probability function p n> k based on (3.1) 
brought no explicit formula for p n ^ for generic k. The only analytic results avail- 
able are due to Edelman (1997) who proved the following properties of the above 
probability function 1 1 : 

Property 1. The probability of having all n eigenvalues real equals 

Pn,n - 2-"<"- 1 )/ 4 . (3.2) 

§ The coefficient S(p, k) is the Stirling number of the second kind. 
|| Table 2 provides a useful illustration of both properties. 
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Property 2. For all < k < n, the p n , k are of the form 

Pn,k = r n ,k + s ntk V2, (3.3) 
where r„ ifc and s n ^ are rational. 

The first result is a simple consequence of (3.1) which, at k = n, reduces to a known 
Selberg integral. All known examples suggest that this is the smallest probability out 
of all p n fc's. The second result is based on more involved considerations of (3.1) and 
(2.4). 

3.2. Joint probability density function of complex eigenvalues ofHG T{n/k) 

To evaluate the sought probability function p n ,k, we start with the definition (3.1). In 
a first step, we carry out the A- integrations therein to assess the j.p.d.f. of all complex 
eigenvalues of a random matrix H. € T(n/fc) having exactly k real eigenvalues: 

Pn k (zi,---,zt) = T\ dX i P- HeT ( n/k) {X 1 ,---,X k ;z l ,---,zt>). (3.4) 

i=1 ^H 

To proceed, we consult (2.5) to spot that a part of it, 

x k , e ({z,z})=f[ f dx t n e- A ? /2 n^-^-i nn^-^-^)' ^ 

i=l jR j=l i>j j=l i=l 

coincides up to a prefactor to be specified below with the average characteristic 
polynomial % 

V k A{z,z})=[ I flOe- tr ° 2 / 2 ' 



r 1 

J 3 = 1 

of a k x fc real symmetric matrix C? = C T drawn from the GOE. More precisely, 

IkA{z, z}) = s k k\ Vk,e({z, z}), 
where s k is given by the Selberg integral (Mehta 2004) 

k k k k 

^=^n/^ri e ^ 2/2 ni A «- A .i= 2fe/2 n r ^/2). 0.7) 

■ i=i jR j=i i>j j=i 

Consequently, the j.p.d.f. of all complex eigenvalues of Ti £ T(n/k) is expressed in 
terms of the average GOE characteristic polynomial V k ,e{{ z , z}) as 

Pn k {zi,---,Zt) = ~lf. —Pn,n VkA{z,z}) TT \Zi - Zj\ 2 \Zi - Zj\ 2 

1 1 - Sn i>j=l 

x n {»i - ~z 3 ) erfc cxp (-^^) • (3-8) 

1 Equations (3.4) and (3.5) suggest that Vo,n = 1- 
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A little bit more spadework is needed to appreciate the beauty hidden in this 
representation. Borrowing the result due to Borodin and Strahov (2005), who 
discovered a Pfaffian formula for a general averaged GOE characteristic polynomial 
(see also an earlier paper by Nagao and Nishigaki (2001)), we may write Vk,e in the 
form 

Sn/Sk 



V ktt ({z,z}) = 



^2t({z,z}) 




(3.9) 



21x21 



Here, V n (zi,Zj) is the so-called D-part of the 2x2 GOE matrix kernel (Tracy and 
Widom 1998) to be defined in (3.12) and (3.13) below; A 2 e({z 7 z}) is the Vandermonde 
determinant + 



A 2 i({z, z}) = A 2 i{z 1 ,z 1 , ■■• ,ze,z e ) = Y\_ \z% - zj\ 2 \z t - zj\ 2 ~[[( z i ~ z i)- ( 3 - 10 ) 



i>j 



i>j 



i=l 



Combining (3.8), (3.9) and (3.10), we obtain: 



P-H k (zi ■ ■ ■ , ze) = 



Pn,r 



x erfc 

3=1 



iV2 



pf 



T^n iz% , Zj ) T^n {z% i Zj ) 
T^n {%% i Zj ) T^n , Zj ) 



(3.11) 



2£x2i 



where (Re zi, ■ ■ ■ , Re^) e M. e and (Im Z\, ■ ■ ■ , Im^) e (K+) £ by derivation. 

Equation (3.11) is the central result of this section. Announced in (1.4), it 
describes the j.p.d.f. of I pairs of complex conjugated eigenvalues {zj,Zj} of an n x n 
random matrix Tik whose k remaining eigenvalues are real, k + 2£ = n. 

To make the expression for the j.p.d.f. Pji k explicit, we have to specify the kernel 
function T> n (zi,Zj). The latter turns out to be sensitive to the parity of n (see, e.g., 
Adler et al (2000)). For n = 2m even, the kernel function is given by 

= I P -(*W)/2 i^-Ay) - i^( x ) i2 j+ i(y) 

h-j 



V 2rn (x,y) = -e 



3=0 



(3.12) 



(3.13) 



while for n = 2m + 1 odd, it equals 

t, ,„ , a _ 1 e -(x 2 + v 2 )/2 v^ 1 g2j+i (g) hi (y) - hi (x) q 2j +i (y) 

^2m+l\X, y) — e ^ , 

3=0 3 

Both representations (3.12) and (3.13) involve the polynomials qj(x) skew orthogonal 
on M with respect to the GOE skew product (Mehta 2004) 

(/, g) = \ f dx e~* 2 / 2 f dy e -« 2 / 2 sgn(y - x) f(x) g(y) (3.14) 



such that 

(?2fe, 92^+1 ) 



(l2k+i,Q2e) — h k S k j, (q2k,Q2e) = (to+i, <?2^+i) — 0. 



(3.15) 



+ Throughout the paper, we adopt the definition Afc(rr) = Yli>j( x i ~ x j) which differs from that of 
Borodin and Strahov (2005) who use the same notation Aj,(ai) for the double product with i < j. 
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The skew orthogonal polynomials Qj(x) can be expressed * in terms of Hermitc 
polynomials as (j 

1 



q 2j {x) 
Q2j+i{x) 



22j 



H 2j {x), 



1 



r [F 2j+1 ( a; )-4i^ 2j _ 1 ( a; )] 



(3.16) 



while "tilded" polynomials q~j(x) f entering (3.13) are related to qj(x) via 

[2j)\ 2 2 " 1 ™ 1 
92," (s) = q2j(x) - ^j2mj\V™( x )> 

q2j+i{x) = q2j+i{x). 
Specifying the normalisation 

hj = {q2j, ?2j+i) = X 
completes our derivation of (3.11). 



22i 



(3.17) 
(3.18) 



3.3. Probability function p n ,k as a Pfaffian integral and inapplicability of the Dyson 
integration theorem for its calculation 

The results obtained in Section 3.2 allow us to express the probability function p n ,k 
in the form of an £-fold integral 

i t 



Pn.k = 



Pn.n 

IT 



d 2 z, 



x erfc 



m Zj >0 

Zi — z 



iV2 



pf 



T^n {%i i Zj ) T^n i Zj ) 
T^n {z~i i Zj ) T) n {Z{ , Zj ) 



(3.19) 



21x21 



involving a Pfaffian. It can also be rewritten as a quaternion determinant 

T> n (z i ,Z 1 ) V n (Zi,Zj) 

V n (z~i,Zj) V n (zi,Zj) 



pf 



qdet[D n (zi, Zj)] e> 



21x21 



where the self-dual quaternion T> n (zi, Zj) has a 2 x 2 matrix representation 
@[V n (zi, zj)] = 



T^n {z% , Zj ) T) n (Zi , Zj ) 
T^n (Zi , Zj ) T^n {Zi , Zj ) 

The Pfaffian/quatcrnion determinant form of the integrand in (3.19), closely 
resembling the structure of both the j.p.d.f. (2.12) of all complex eigenvalues and 
the p-point correlation function (2.16) in GinSE, makes it tempting to attack the £- 
fold integral (3.19) with the help of the Dyson integration theorem (see Section 2.3). 
Unfortunately, the key condition of this theorem - the projection property - is not 
fulfilled. 

To see this point, we represent the kernel function T> n (x,y) in the form 

-. 71— 1 

V n {x,y) = - e-(* 2 +« 2 )/ 2 qjWfijkQkiv), (3-20) 

j,k=0 

* The representation (3.16) is not unique; see, e.g., Eynard (2001). 
(J Equation (3.16) assumes that H-\(x) = 0. 

f Note that the q2j(x) is no longer a polynomial of the degree 2j. 



3 Statistics of real eigenvalues in GinOE spectra 



22 



see the primary definition (3.12) and (3.13). In (3.20), the real antisymmetric matrix 
fi of the size n x n depends on the parity of n. It equals 



and 



\ 



(3.21) 



-odd 



^ m — 1 



(3.22) 



\ £o ••• / 

for n = 2m and n = 2m + 1, respectively. We remind that hj is defined by (3.18); also 
we have used the notation 



e 2 



1 

-1 



1 



m! 



(3.23) 



Actually, the representation (3.20) can be put into a more general form due to Tracy 
and Widom (1998) that would contain arbitrary, not necessarily skew orthogonal, 
polynomials upon a proper redefinition of the matrix /}. 

For the Dyson integration theorem to be applicable, the projection property for 
the self-dual quaternion T> n (zi,Zj) must hold. For this to be the case, the integral 
identity 



/ 



da(w) [V n (z 1 ,w)T> n (w,z 2 ) - T> n (z 2 ,w)T> n (w, zi)] = -T> n (z 1 ,z 2 ) 



should be satisfied for the measure 

' w — w 



da{w) = erfc 



iV2 



9(lmw) d 2 w. 



(3.24) 
(3.25) 



Here, 9((f>) is the Heaviside step function, 

1 if <p > 



6{4>) 



if <b< 



Straightforward calculations based on (3.20) show that the integral on the l.h.s. of 
(3.24) equals 



± e -(*?+*2)/2 

2 



1 



j,k=0 



(3.26) 



where n x n matrix % has the entries 

Xjk = \j da(w)e-^ 2+ ^ 2 [q,(w)q k (w) - q 3 (w)q k (w)] . (3.27) 

Since | (/x%) + -in, the l.h.s. of (3.24) given by (3.26) differs from the r.h.s. (3.24). 
As a result, the kernel function I?„(zi,Z2) does not satisfy the projection property 
§. Consequently, the Dyson integration theorem is inapplicable for the calculation of 
■pnM in the form of the Pfaffian integral (3.19). 

% That —(Ax) cannot be a unit matrix 1, follows from the fact that \jk ls purely imaginary [(3.27)] 
while fijk are real valued [(3.21) and (3.22)]. See Appendix B for an explicit calculation. 
§ As soon as the integral on the l.h.s. of (3.24) combines a D-part of the GOE 2x2 matrix kernel 
originally introduced for the GOE's real spectrum with the GinOE-induccd measure da(w) supported 
in the complex half-plane lm.w > 0, a violation of the projection property is not unexpected. 
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4. Probability function p n ,k'- Sensing a structure through particular cases 

Before turning to the derivation of the general formula for p n ^ (see the results 
announced in Section 1.2), it is instructive to consider a few particular cases 
corresponding to low values of £, the number of pairs of complex conjugated 
eigenvalues. Below, the cases of I = 1, 2 and 3 are treated explicitly. 



4-1. What is the probability to find exactly one pair of complex conjugated 
eigenvalues? 

As a first nontrivial application of the Pfaffian integral representation (3.19) for the 
probability function p n ,k, let us consider the next-to-the-simplest || case of £ = 1 
corresponding to the occurrence of exactly one pair of complex conjugated eigenvalues. 
Since the £)-kernel (3.20) is antisymmetric under exchange of its arguments, 

V n (x,y) = -V n (y,x), (4.1) 

the Pfaffian in (3.19) reduces to 

V n (z,z) ' 

-V n {z,z) 



pf 

resulting in 



T> n (z,z) 



Pn,n-2 = Pn,n T / d 2 Z Crfc ( — =f ) T> n (z, z). (4.2) 



lmz>0 



iV2 



To calculate the integral 



Ji= / d 2 z erfc ( ) V n (z,z) = da(z)V n (z,z) 

Jlmz>0 V W2 / J 

(see (3.25)), we rewrite it in a more symmetric manner 
i J da(z) V n {z,z) -V n {z,z) , 

and make use of (3.27), (3.25) and (3.20) to deduce that it equals 

1 ™ _1 1 
~9 VjkXk] =-^tr (0 ,„_i)(Ax), 



j,k=0 

or, cquivalently, 



2i 



J da(z) T> n (z, z) = j tr (0in _i)CT. (4.3) 

Here, a is given by a = lifix ( see a ls° Appendix B). We therefore conclude that the 
probability sought equals 

Pn,n-2 = ^Pn,ntr( 0)Tl _i)d\ (4.4) 

Due to the trace identity (C.5) proven in Appendix C, we eventually derive: 

Pn,n-2 = Pn.n tr (0) [n/2J -1) (4-5) 

reducing the size of the matrix by two. Here, the smaller matrix q depends on the 
parity of n, as defined in the Section 1.2 (see also (C.3) and (C.4)). 

|| In the simplest case of I = 0, our representation (3.19) reproduces the result (3.2) first derived by 
Edelman (1997). 
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Remarkably, the trace in (4.5) can explicitly be calculated (see Appendix D) 
to yield a closed expression for the probability to find exactly one pair of complex 
conjugated eigenvalues: 



Pn,n-2 = 2p n ,n / dy y e y2 erfc(y\/2) L 2 _ 2 (-2y 2 ). 
Jo 



(4.6) 



Yet another, though equivalent, representation for the probability p n , n -2 is given in 
Section 7 that addresses the large-n behaviour of p n ,n-2- 



4-2. Two pairs of complex conjugated eigenvalues (1 = 2) 
For 1 = 2, the Pfaffian in (3.19) 

V n {zi,z{) V n (z 1 ,z 2 ) T> n {zi,z 2 ) 



pf 



-T> n (zi,zi) V n (zi,z 2 ) V n (zi,z 2 ) 

-T> n (z 1 ,z 2 ) -T> n (z- L ,z 2 ) V n (z 2 ,z 2 ) 

-V n (z l7 z 2 ) -V n (z l7 z 2 ) -V n {z 2 ,z 2 ) 

reduces to 

T> n (zi,zi)V n (z 2 ,z 2 ) +T> n (zi,z 2 )V n (zi,z 2 ) - V n {z\,z 2 )V n {z\,z 2 ), 
so that 

Pn,n-4 = -2 p n ,n Jda(z!) J dct(z 2 ) 

x V n (zi,zi)V n (z 2 ,z 2 ) +V n (zi,z 2 )V n (zi,z 2 ) - V n {zi,z 2 )V n {zi,z 2 ) 
Apart from a known integral taking the form of (4.3), a new integral 
J 2 = da(zi) da(z 2 ) T> n (zi, z 2 ) V n {z\, z 2 ) - V n {z\, z 2 ) T> n (zi, z 2 ) 



(4.7) 



(4.8) 



(4.9) 



(4.10) 

appears in (4.9). Somewhat lengthy but straightforward calculations based on (3.27), 
(3.25) and (3.20) result in 



I2 = -tr^^^o- 2 ). 
Combining (4.9), (4.3), (4.10) and (4.11), we obtain: 



Vn.n—A: Pji,n 



1 I 1 I ~2\ 

g (t^Cn-i) 0- ) - ^ tr (o,n-i)(cr ) 



Due to the trace identity (C.5) proven in Appendix C, the latter reduces to 



_ 1 



(tr( ,|„/2J-l)£) - tr( ,Ln/2J-l)(e ) 



(4.11) 
(4.12) 
(4.13) 



Interestingly, the expression in the parenthesis of (4.13) coincides with Z( 1 2)(p 1 ,p 2 ) 
after the substitution tr g 1 = pj (see Table 1). 

4-3. Three pairs of complex conjugated eigenvalues (i = 3) 

The complexity of the integrand in (3.19) grows rapidly with increasing I. For £ = 3, 
that is three pairs of complex conjugated eigenvalues in the matrix spectrum, the 
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T> n (z 1 ,z 3 ) T> n (zi,z 3 ) \ 

1>n(zi,z 3 ) T> n (zi,z 3 ) 

T^n(Z2,Z 3 ) T> n (z 2 ,Z 3 ) 

T> n {z 2 ,Z 3 ) V n (z 2 ,Z 3 ) 

V n (z 3 ,z 3 ) 
-V n {z 3 ,z 3 ) J 



Pfafhan of the 6x6 antisymmetric matrix 

/ T> n {zi,z{) V n (z 1 ,z 2 ) V n {z\,z 2 ) 

-V n (z\,z\) X>„(zi,z 2 ) T> n (zi,z 2 ) 

-V n (z-i,z 2 ) -T> n (zi,z 2 ) T) n (z 2 ,z 2 ) 

-V n (zi,z 2 ) -V n (z 1} z 2 ) -V n (z 2 ,z 2 ) 

-V n {zi,z 3 ) -T> n (zi,z 3 ) -V n (z 2l z 3 ) -T> n (z 2 ,z 3 ) 

\ -T> n (zi,z 3 ) -T> n (zi,z 3 ) -T> n (z 2 ,z 3 ) -T> n (z 2 ,z 3 ) 

is getting involved. It can be calculated with some effort to give 15 terms which can 
be attributed to three different groups. The first group consists of the single term 

d =V n (z 1 ,z 1 )V n {z 2 ,z 2 )V n {z 3 ,z 3 ). (4.14) 

The second group contains 6 terms, 

G 2 =T> n (z 1 ,z- L ) V n (z 2 ,z 3 )V n (z 2 ,z 3 ) -V n {z 2 ,z 3 )T> n {z 2 ,z 3 ) 
+ V n (z 2 ,z 2 ) r D n (z 1 ,z 3 )T> n (z 1 ,z 3 ) -V n (z 1 ,z 3 )V n (z 1 ,z 3 ) 

+ T> n (z 3 ,z 3 ) T> n (zi,z 2 )V n (zi,z 2 ) -V n (zi,z 2 )V n (zi,z 2 ) , (4-15) 
while the third group counts 8 terms: 

G 3 = T> n (zi,z 3 ) T> n (z 2 , z 3 ) T> n (zi, z 2 ) - T> n (zi, z 3 ) T> n (z 2 , z 3 ) T> n (zi,z 2 ) 
+ T> n (zi,z 2 ) V n (z 2l z 3 ) T> n (zi,z 3 ) - T> n (zi, z 2 ) V n (z 2l z 3 ) £>„(zi, z 3 ) 
+ V n (zi,z 2 ) V n (z 2 ,z 3 ) V n (z 1} z 3 ) - V n (z-i,z 2 ) T> n (z 2 ,z 3 ) V n (z 1} z 3 ) 
+ V n (zi,z 2 ) V n (z 2 ,z 3 ) T> n {zx,z 3 ) - X>„(zi, z 2 ) T> n (z 2 ,z 3 ) V n (z 1 ,z 3l ). 
In the above notation, the probability function p n , n -6 takes the form 
4i 

Pn,n—6 

The integrals containing Gi and G 2 can easily be performed with the help of (4.3), 
(4.10) and (4.11) to bring 

3 

(4.18) 



■Pn,n J da{zi) J da(z 2 ) J da(z 3 ) G x +G 2 + G 3 



(4.16) 
(4.17) 



J ' da(zi) J da(z 2 ) J da(z 3 ) Gi = T\ = ^ tr( , n _i)d- 
and 

J da(zi) j da(z 2 ) J da(z 3 ) G 2 = 3 X\ X 2 



3 -tr 



(0,n-l) ( 



1 

-tr 



(0,n-l 



)(^ 2 ) 



(4.19) 



The remaining integral involving G 3 can be evaluated similarly to T\ and I 2 , the result 
being 

l 3 =Jda(z 1 ) J da(z 2 ) J da(z 3 ) G 3 = ^ tr( 0i „_i) (<x 3 ). (4.20) 
Combining (4.17), (4.18), (4.19) and (4.20), we derive: 

Pn,n-6=Pn,n ^ (t r (0,n-l) CT J ~ g t r (0,n-l)°" tr (0,n-l) ( CT ) + g tr (0,n-l) ( CT ) • 

(4.21) 
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Finally, we apply the trace identity (C.5) to end up with the formula 



Pn.n — 6 — „ Pn,n 

6 



(tr(o,|n/2j-i)f?) 



- 3tr (0! |„/2j_i)f?tr(o,|„/2j-i)(£ 2 ) + 2tr (0 ,Ln/2j-i)(£ 3 ) . (4.22) 

The expression in the parenthesis of (4.22) is seen to coincide with -Z(i3)(pi,p2,.P3) 
after the substitution tr q 3 = pj (see Table 1). 

4-4- Higher I 

The three examples considered clearly demonstrate that the calculational complexity 
grows enormously with increasing £, the number of complex conjugated eigenvalues in 
the random matrix spectrum. Indeed, the number of terms Afe in the expansion of the 
Pfafhan in (3.19) equals Mi = (21 — 1)!! and exhibits a faster-than-exponential growth, 
Afe ~ 2 l+1 / 2 e ^( ln ^- 1 ) j for £ ^> 1. For this reason, one has to invent a classification of 
the terms arising in the Pfaffian expansion that facilitates their effective computation. 
This will be done in Section 5, where we introduce a topological interpretation of the 
Pfaffian expansion, and prove the Pfafnan integration theorem which can be viewed 
as a generalisation of the Dyson integration theorem (see Section 2.3). 



5. Topological interpretation of the Pfaffian expansion, and the Pfaffian 
integration theorem 

5.1. Statement of the main result and its discussion 

Before stating the main result of this section, the Pfaffian integration theorem, we 
wish to start with presenting a simple Corollary to the Dyson integration theorem. 

Corollary 5.1. Let f{x,y) be a function with real, complex or quaternion values 
satisfying the conditions (2.8a) and (2.8b) of the Theorem 2.1, and dir be a suitable 
measure. Then 

J f[ dn( Xj )det [f{x l ,x j )} £xe = r ^|^ » ( 5 - la ) 

where 

c = y dn(x) f(x, x). (5.1o) 

For f taking quaternion values, the det should be interpreted as qdet, the quaternion 
determinant (Dyson 1972). 

Proof. Repeatedly apply the Dyson integration theorem to the l.h.s. of (5.1a) to 
arrive at its r.h.s. ■ 



Importantly, the above Corollary exclusively applies to functions f(x,y) satisfying 
the projection property as defined in Section 2.3 (see Definition 2.1 therein). How- 
ever, guided by our study of the integrable structure of GinOE, we are going to ask if 
the integrals of the kind (5.1a) can explicitly be calculated if the projection property is 
relaxed. In general, the answer is positive. In particular, for f(x,y) being a self-dual 
quaternion, the following integration theorem will be proven. 
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Theorem 5.1 (Pfaffian integration theorem). Let dn(z) be any benign measure 
on z <G C, and the function Q n {x, y) be an antisymmetric function of the form 

_^ n—l 

Qn(x,y) = - ^2 qj(x) fijkqkiy) (5.2a) 



j,k=o 



where the qj(x) are arbitrary polynomials of j-th order, and fl is an antisymmetric 
matrix. Then the integration formula 



n dn ( z ^ pf 

i=i 



Qn {z% , Zj ) Q n (Zi , Zj ) 

Qn (Zi , Zj ) Qn (Zi , Zj ) 

i\ £ ( \ 

2) V) ( 9 tr (0,n-l)*>" 



2£x2^ 



1 

-tn, 



.„ „ ■ . , 2 «r(o,n-i)V J ( 5 - 2& ) 

holds, provided the integrals in its l.h.s. exist. Here, are zonal polynomials whose 

£ arguments are determined by a matrix v with the entries 

n—l 

v a ,fi = Aa.fe / d7T(» [^(z) ^(z) - g^z) . (5.2c) 



fe=0 



As we integrate over all variables, the Pfaffian integration theorem can be viewed 
as a generalisation of the Corollary 5.1 proven a few lines above, for the case of a 
kernel not satisfying the projection property. This follows from the identity 

Qn{zi, Zj) Qn(zi, Zj) 
Qn{zi, Zj} Qn{zi, Zj) 



Pf 



= qdet 



2£x2£ 



Qn{Zi, Zj) 



where the quaternion Q n (zi, Zj) has the 2x2 matrix representation: 



f(zi,Zj) = Q[Q n (zi,Zj 



QrJ\Zi 7 Zj ) 
Qn {%i ■> Zj) 



Qn 1 Zj ) 
Qn {Zi, Zj ) 



(5.3) 



To see that the Pfaffian integration theorem reduces to the Corollary for the particular 
case of a kernel with restored projection property, we spot that the latter is equivalent 
to the statement 



v = -2i 1„ 

as can be deduced from the discussion below (3.24), Section 3.3. As a result, 

/2 X 3 
tr( ,„-i)t> J =n I 

and the r.h.s. of (5.26) reduces to (Macdonald 1998) 

-i(-2i*r 



n ) V) It 



(1 + z)"/ 2 
Finally, noticing from (5.3) that 




z=0 



J dir(z)f(z,z) = J dn(z) 



(z,z) 

Q n (z,z) 



(5.4) 
(5.5) 



z=0 



(5.6) 



■ eotr ( ,„-i)t> = (n/2) e 



(5.7) 
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Table 3. The vocabulary of topological terms defined to interpret the ordered 
Pfaffian expansion. The notation used is: D — Definition, E - Example, L - 
Lemma, T — Theorem, C - Corollary, F - Figure. 


Term 


Notation 


Appearance 


String 


s t 


D5.1, E5.1, F3, F4 


Length of a string 


IIS, || 


D5.2, E5.2 


Equivalent strings 


Si ~ Sj 


D5.3, E5.3, L5.1, F3 


Equivalence class of strings 




D5.4, E5.4, L5.1, F3, F5 


Size of equivalence class 


lie, II 


D5.4, L5.1, F3 


Substring 




D5.5, E5.5 


Length of a substring 


n - ||c(p)|| 




Loop-like substring 


s<*> 


D5.6, E5.6, L5.2, L5.5, 
F3, F4 


Longest loop-like string 




F3, E5.7, L5.4 


Adjacent loop-like string 


S« 


D5.7, E5.7, E5.9, E5.ll, 
L5.4-L5.6, C5.2, F4, F5 


Handedness of adjacent (sub)string 


M(a L ,ce R ) 


D5.8, E5.8, E5.9, E5.ll, 
L5.6, C5.2, F5 


Equivalent adjacent (sub)strings 


Si ~ Sj 


D5.9, E5.10, E5.ll, 
L5.7, F5 


Equivalence class of adjacent (sub)strings 


ACj 


D5.10 


Compound string 


s, 


D5.ll, T5.3, F3, F6 


Topology class 


A 


F6 



with e = diag(l, 1), we conclude that the constant c in the Corollary [(5.16)] equals 
c = n/2 so that the result (5.6) brought by the Pfaffian integration theorem is 
identically equivalent to the one [(5.1a)] following from the Dyson integration theorem. 
We stress that this is only true for the kernel Q n (zi,Zj) satisfying the projection 
property. 

To prove the Theorem 5.1, we will invent a formalism based on a topological 
interpretation of the ordered Pfaffian expansion. For the readers' benefit, a vocabulary 
of the topological terms to be defined and used in the following sections is summarised 
in Table 3. 

5. 2. Topological interpretation of the ordered Pfaffian expansion 

To integrate the Pfaffian in (5.26), we start with its ordered expansion 



pf 



Qn {%i i %j ) Qn i %j ) 



Qji i Zj ) Qn {%% i %j ) 



= WT\ H Sgn ( CT ) \\Qn{Wa(2j+l),W a{2j+2) ). 

2£x2£ z l - aeS2e j=Q 



(5.8) 

Here, the summation extends over all permutations a E S 2 i of 2£ objects 

{wi = zi,w 2 = zi,- ■ -,W2t-i = z e ,w 2 e = z e } (5.9) 
so that the total number of terms in (5.8) is {2t)\. 
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5.2.1. Strings and their equivalence classes 

Definition 5.1. Each term of the ordered Pfaffian expansion is called a string. 
The i-th string Si equals 

e-i 

Si = SgnOj) Y[ Qn{Wa i (2j+l),W<r i (2j+2))- (5.10) 
3=0 

where a is the i-th permutation out of {2tj\ possible permutations a G Sit- Notice 
that a sign is attached to each string. 

Example 5.1. The ordered expansion of the Pfaffian for I = 2 [see (4.7)] contains 
4! = 24 strings that we assign to three different groups (their meaning will become 
clear below): 









+ (11)(22) 


+(12)(12) 


-(12)(12) 


-(11)(22) 


-(12)(21) 


+(12)(21) 


+(11)(22) 


+(21)(21) 


-(21)(21) 


-(11)(22) 


-(21)(12) 


+(21)(12) 


+(22)(11) 


+(12)(12) 


-(12)(12) 


-(22)(11) 


-(12)(21) 


+(12)(21) 


+(22)(11) 


+(21)(21) 


-(21)(21) 


-(22)(11) 


-(21)(12) 


+(2T)(12) 



For brevity, the obvious notation ±(pq)(pq) was used to denote the string 
±Qn(z p , z q )Qn(z p , z q ). The three strings shown in bold are those that previously 
appeared in (4.8) when treating the probability p n , n -i- 

Definition 5.2. The length ||<Sj|| of a string Si equals the number of kernels it 
is composed of. 

Example 5.2. The string Si in (5.10) is of the length I: \S { \ = I. All strings 
in (5.11) are of the length 2. 

Definition 5.3. Two strings Si and Sj of the ordered Pfaffian expansion are said 
to be equivalent strings, Si ~ Sj, if they can be obtained from each other by (i) 
permutation of kernels and/or (ii) permutation of arguments inside kernels (these will 
also be called intra-kernel permutations). 

Example 5.3. For 1 = 2, three different groups of equivalent strings can be identified 
as suggested by (5.11). The first group, exemplified by the string +(11)(22) involves 
8 equivalent strings; two other groups, each consisting of 8 strings as well, are repre- 
sented by the strings +(12)(12) and — (12)(12), respectively. 

Definition 5.4. A group of equivalent strings is called the equivalence class of 

strings. The j-th equivalence class to be denoted as Cj consists of \\Cj\\ equivalent 
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Ordered Pfaffian Expansion: 

[21 — 1)!! equivalence classes 
each made of (2£)\\ equivalent 
strings 




Longest Loop-Like Strings 

Total amount: (2£)V.(2£ - 2)!! 



Compound Strings 

made of at least two 
sets of loop-like strings 



Figure 3. The (2£)\ terms, or strings, of an ordered PfafBan expansion can be 
assigned to (21 — 1)!! equivalence classes {Ci, • • • ,C(2£— i)!l}> each one containing 
(2£)\\ equivalent strings (Lemma 5.1). Any given string of the length I can be 
decomposed into a set of loop-like substrings of respective lengths {£j } such that 
2~2j = ^ (Lemma 5.2). The sets consisting of only one loop-like string are called 
longest loop-like strings; the sets made of more than one loop-like substring are 
called compound strings. An amount of longest loop-like strings is counted in 
Lemma 5.3. 



strings. The number \\Cj\\ is called the size of equivalence class. 

Example 5.4. For I = 2, exactly three different equivalence classes of strings (5.11) 
can be identified in the ordered Pfaffian expansion (5.8). 

In the context of the initiated classification of strings arising in the ordered Pfaffian 
expansion (5.8), a natural question to ask is this: Can the total number of equivalence 
classes and the number of of equivalent strings in each class be determined? The 
answer is provided by Lemma 5.1. 

Lemma 5.1. All terms of the ordered Pfaffian expansion can be assigned to (21 — 1)!! 
equivalence classes {C\, ■ ■ ■ ,C^e-iy.\}, eac h containing (2i)\\ equivalent strings: \\Cj\\ = 
(2£)\\ for all j e 1, • • • , (21 - 1)!!. 

Proof. Consider a string Si belonging to the equivalence class Cj and composed 
of £ specific kernels, ||<Sj|| = I. There exist £\ possible permutations of kernels and 
2 l intra-kernel permutations of arguments. As a result, the total number of strings 
generated by these two operations from a given string Si equals ||Cj|| = 2 l £\ — (2£)\\. 
Since the total number of strings in the ordered expansion is (2£)\, one concludes that 
there exist (2£)\/(2£)\\ = (21 — 1)!! different equivalence classes. ■ 

The results of this and the two subsequent Sections are summarised in Fig. 3. 
5.2.2. Decomposing strings into loop-like substrings 

Having assigned all (2£)\ strings of the ordered Pfaffian expansion to (2£— 1)!! equiv- 
alence classes each containing (2^)!! strings, we wish to concentrate on the structure 
of the strings themselves. Below, we shall prove that any string Si (of the length 
||<Si|| = £) can be decomposed into a certain number (between 1 and £) of loop-like 
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substrings, see the Lemma 5.2. To prepare the reader to the definition of a loop-like 
substring, we first define the notion of a substring itself. 

Definition 5.5. A product of p kernels Q n is called a substring of the string 

e-i 

S i = sga(a i )Y[Q n (w (Ti (2 j+1 - ) ,w <Ti (2j+2)), W^iW = t, ( 5 - 12 ) 

3=0 

specified in Definition 5.1, if it takes the form 
v 

S i P) = II Qn(w^(2 3k +l),W^ (2jk+2) ) (5.13) 
fe=l 

where ji j 2 ^ ■•• ^ j P - The length \\S^\\ of the substring is \\S^\\ = p with 
1 < p < I. Notice that no sign is assigned to a substring. 

Example 5.5. The string 

+(11)(23)(23) 

arising in the I = 3 Pfaffian expansion (see the first term in (4.15)) can exhaustively 
be decomposed into seven substrings % 

(11), (23), (23), (11)(23), (11)(23), (23)(23), (11)(23)(23) 

V v ' V v ' V v ' 

p=l p=2 p=3 

of lengths p = 1, 2 and 3, respectively. 

Definition 5.6. A substring 

v 

$i P) = Y[Qn(w ai{2 j k + l),W ai{2jk+ 2)) 
fe=l 

is said to be a loop-like substring of the length p = ||«S t ^||, if the two conditions are 
satisfied: 

(i) The set of all 2p arguments 

W 2p = {w CT .( 2 j 1+ l),W (7 .(2j 1+ 2), • • • ,W a .(2j p+ l), W a .(2j p +2)} 

collected from the substring s\ p ^ remains unchanged under the operation of 
complex conjugation 

yV 2p = {w ad2n + 1) , W ad2n+2) , ■ ■ ■ ,«V 4 (2j p +l),«V 4 (2j p +2)}- 

That means the two sets W 2p and W 2p of arguments are identical, up to their 
order. (This property will be referred to as invariance under complex conjugation.) 

(ii) For all subsets 6S^ consisting of q kernels Q n with 1 < q < p— 1, the substring 
S^ P \8S^ of the length p — q obtained by removal of 5S^ from is not in- 
variant under the operation of complex conjugation of its arguments. 



1 It is easy to see that the number of substrings of the length p equals ( ) so that the total amount 
of all possible substrings of a string of the length £ is 



T- - 1. 
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Example 5.6. Out of seven substrings of the string +(11)(23)(23) detailed in the 
Example 5.5, the two substrings 

(11), (23)(23) 

are loop-like (of the lengths p = 1 and 2, respectively). The remaining five substrings 
are not loop-like. Four of them, 

(23), (23), (H)(23), (11)(23), 

are not loop-like substrings because the property (i) of Definition 5.6 is not satisfied. 
The fifth substring (of the length p = 3) 

(1T)(23)(23) 

is not loop-like because the property (ii) of Definition 5.6 is violated. Indeed, there 
does exist a subset SS^ consisting of one kernel, represented by the pair of argu- 
ments (11), whose removal would not destroy the property (i) for the reduced substring 
(23)(23). 

The example presented shows that a particular string of the ordered Pfaffian expan- 
sion could be decomposed into a set of loop-like substrings. Is such a decomposition 
possible in general? The answer is given by the following Lemma. 

Lemma 5.2. Any given string Si of the length \\Si\\ — I from the ordered Pfaf- 
fian expansion can be decomposed into a set of loop-like substrings Sf 3 ^ of respective 
lengths =lj, 

3 

such that J2j £j — 

Proof. We use induction to prove the above statement. 

(i) Induction Basis. For £ — 1, the Lemma obviously holds since the strings (1,1) 
and (1,1) are loop-like by Definition 5.6. 

(ii) Induction Hypothesis. The Lemma is supposed to hold for any string Si of the 
length \\Si || =1: 

Si = \JS^\ with ^2lj=l (5.14) 

3 3 

(iii) Induction Step. Consider a given string Si of the length ||<Sj|| = I + 1. Given the 
induction hypothesis, we are going to prove that such a string Si can be decom- 
posed into a set of loop- like substrings. 

To proceed, we note that any given string Si of the length ||<Sj|| = I + 1 can 
be generated from some string Si of length I (see 5.12) by adding to it an addi- 
tional pair of arguments (z£+i,Z£+i): 

£ 

Si = (Z£ +1 ,Z£ +1 ) ® J|^( 2 ^+ 1 )' W ^( 2 ^+ 2 )) ( 5 ' 15 ) 

fc=l 

S v ' 

the string Si with sgn(cri) dropped 

with (or without) further exchange of either zp,+i or zi + \ with one of the arguments 
belonging to the string Si of length I. 
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(a) If no exchange is made, the given string <Sj is a unit of a single loop-like string 
n^ 1 -* = (zt+i,zt+i) and of a string Si admitting the decomposition (5.14). 
As a result, the string Si of the length (£ + 1) is decomposed into a set of 
loop-like substrings 

S^U^U^jS^y (5.16) 

(b) If either zi + \ or zg + i was swapped with one of the arguments belonging to 

(£ ■ ) 

a loop- like substring S i 3 " C Si of the string Si , such an exchange will give 

(£'■ ) 

rise to a new loop-like substring S i 30 C Si of the length 

ll^ 0) || = ll^ll + 1 = ^ + 1. 

Consequently, the given string Si of the length (£ + 1) is then decomposed 
into a set of loop-like substrings 

(£'. ) 

To prove that the substring S { 30 is indeed loop-like, two properties have to 
be checked in accordance with the Definition 5.6. 

First, one has to show that the set of 2(£j + 1) arguments collected from the 

(p. ) 

substring S { 30 is invariant under the operation of complex conjugation; this 

(£ ■ ) 

is obviously true because S t 30 is loop-like. Second, one has to demonstrate 
that the removal of any subset 5S from S^ will destroy the invariance 

(£' ) 

property of the remaining substring S i 30 \SS. Three different cases are to 
be considered here: 

(bl) If the subset SS does not contain the fragments (2^+1, •••) and 

(•••,Z£ + i), it is also a subset of Si . Since the latter is loop-like, the 
invariance property is destroyed. 

(b2) If the subset 5S contains only one of the fragments (2^+1, • • •) or 
(• • • , zi+i), the invariance property is obviously destroyed. 

(b3) If the subset SS contains both fragments (zt+i,u) and (v,zi + i), the 
invariance property is also destroyed. To prove it, we use the reductio ad 

(£'. ) 

absurdum. Indeed, let us assume that there exists a subset 5S C S i 30 C Si, 
containing both fragments (zt+i,u) and (v,zt+i), whose removal does not 

destroy the invariance property of S i 30 \SS. The existence of such a subset 

(£ ■ ) 

SS implies the existence of yet another subset SS' C Si 30 C Si, 

SS' = {SS\ {(z i+1 ,u) U (v, z t+1 )} } U {(v, «)}, (5.18) 

(£■ ) 

whose removal does destroy the invariance of S t 30 \SS' under the operation 
of complex conjugation (this claim is obviously true because the substring 
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Figure 4. Graphic explanation of the notion "loop-like substring", (a) The 
substring (23) (23) can be transformed into the form of an adjacent substring 
(in the sense of Definition 5.7) by flipping arguments in the second kernel, 
(23) i— > (32). (b) The emerging adjacent substring (23)(32) can clearly be depicted 
in the form of a loop by gluing the arguments 2 and 2 together. 



S\ i0> is loop-like). Then, the identity 

S { f io) \8S = sf io) \5S l (5.19) 

suggests that the removal of 5S from s\ J0 ' must destroy the invariance 

(p. ) 

of 5^ 30 \5S as well. Since this contradicts to the assumption made, one 
concludes that the invariance property is indeed destroyed. 

End of the proof. ■ 



We close this Section by providing a brief explanation of the origin of the term "loop- 
like substring" introduced in Definition 5.6. In Section 5.2.4, it will be proven that any 
loop-like substring can be brought to the form of an adjacent substring (Definition 5.7) 
by means of proper (i) permutation of kernels and/or (ii) permutation of intra- kernel 
arguments (Lemma 5.5). The latter can graphically be represented as a loop. Indeed, 
the loop-like substring (23) (23) considered in Example 5.6 can be transformed into the 
form of an adjacent substring (in the sense of Definition 5.7) by flipping arguments in 
the second kernel, (23) i— > (32): 

(23) (23) i — ► (23) (32) 

loop — like substring adjacent substring 

The resulting adjacent string (23) (32) can be drawn in the form of a loop, see Fig. 4. 
This is precisely the reason why the substring (23) (23) is called loop-like. 

5.2.3. Counting longest loop-like substrings S\ of the length i 

Although in this subsection, we are going to concentrate on the longest loop-like 
substrings + of the length £, our main counting result given by Lemma 5.3 stays valid 
for loop-like substrings of a smaller length 1 < p < I. 

Lemma 5.3. The ordered Pfaffian expansion contains (2£)\\ (21 — 2)!! longest loop-like 
strings of the length L 

+ As soon as substrings of the longest possible length I are considered, they are strings themselves. 
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Proof. Let Nl{£) be the total number of loop-like strings of the length £ in the 
ordered Pfaffian expansion and let Ul{£) denote the number of equivalence classes all 
longest loop-like strings can be assigned to. Following Lemma 5.1, the two are related 
to each other as 

M L (£) = (2£)\\n L (£), (5.20) 

because each equivalence class Cj contains precisely (2£)\\ equivalent strings (see 
Definition 5.3). It thus remains to determine nj,(£) that can equally be interpreted as 
a number of inequivalent * longest loop-like strings of the length £. 

The latter can be evaluated by counting the number of ways, n^(£ + 1), all 
inequivalent longest loop-like strings of the length £ + 1 can be generated from those 
longest of the length £. Since both numbers, Ul{£) and ul{£ + 1), refer to the longest 
loop-like strings, the two correspond to the Pfaffians of matrices of the size £ x £ and 
(^+1) x (£+1), correspondingly. In the language of strings, an increase of the matrix 
size by one leads to the appearance of an additional pair of arguments (ze+i,ze+i) in 
a string of the length £ + 1 . 

We claim that 

n L (£+l) = 2£n L (£). (5.21) 

To prove this, we concentrate on a given longest loop-like string of the length £ and 
add to it an additional pair of arguments (z^+i, ze+i): 

e 

(z e+ i,z e+ i) O Y[(w az{2jk+1) ,w ad2jk+2) ) . (5.22) 
fe=i 

V v ' 

longest loop — like string of the length £ 

Here, Oi is a particular permutation of 2£ arguments (5.9) corresponding to a longest 
loop-like string of the length £. The resulting string (5.22) is not a longest loop-like 
string of the length £+1 (rather, it is composed of two loop-like strings of the lengths 
1 and £, respectively). Since a loop-like string necessarily assumes the presence of 
a fragment (ze+i, ■■■)("■■> ze+i) somewhere in the string, one has to exchange either 
Z£ + i or Z£ + i with one of the arguments belonging to the original longest loop-like 
string of the length £. Clearly, there exist 2£ exchange options for each argument, Z£+i 
(or ze+i). As a result, we arrive at the relation (5.21). Given m,(l) = 1, we derive 
the desired result by induction: 

n L {£) = {2£ - 2)\\ (5.23) 

Combining it with (5.20) completes the proof, jj ■ 

5.2.4- Adjacent vs non-adjacent loop-like substrings 

Further classification of loop-like substrings is needed in order to prepare ourselves to 

* In view of Definition 5.3, the two strings are inequivalent if they cannot be reduced to each other 
by means of (i) permutation of kernels and/or (ii) intra- kernel permutation of arguments. 
J It is instructive to turn to Example 5.1 that discusses the ordered Pfaffian expansion for i = 2. The 
Lemma 5.3 predicts existence of A/j,(2) = 4!! 2!! = 16 longest loop-like strings that can be assigned to 
n^(2) = 2!! = 2 equivalence classes, each composed of (2£)\\ = 4!! = 8 equivalent strings. This is in 
line with direct counting (5.11): the two equivalence classes are represented by the longest loop-like 
strings +(12)(12) and — (12)(12) (see the second and third column); each equivalence class contains 
8 equivalent strings. 
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the proof of the Pfaffian integration theorem. 

Definition 5.7. A loop-like substring S- p ^ of the length p = \\S^\\ is called ad- 
jacent loop-like substring, or simply a loop, if it is represented by a product 

v 

S ( i P) = II ( 9™( U7 ^(2^+l) J W CTl (2 Jfc + 2)) 
fe=l 

of p kernels such that: 

(i) The first argument of the first kernel and the second argument of the last, p-th 
kernel, are complex conjugate of each other, 

UV i( 2j 1+ l) = W a . {2jp+ 2)- 

(ii) For each pair of neighbouring kernels in the string, the second argument of the 
left kernel in the pair and the first argument of the right kernel in the pair are 
complex conjugate of each other, 

w ai (2j k +2) = w a . { 2j k+1 +i), k = !,■■■, V- 1- 



Example 5.7. Out of 16 longest loop-like strings arising in the Pfaffian expansion 
for I = 2 f, the following eight are adjacent: 

-(12)(21), -(21)(12), -(12)(21), -(21)(12), 
+ (12)(21), +(21)(12), +(12)(21), +(21)(12). 

Notice that although longest loop-like strings have been considered in the above exam- 
ple, the notion of an adjacent string is equally relevant for a loop-like string of smaller 
length. 

Lemma 5.4. Out of (2£)l\(2£ — 2)!! longest loop-like strings of the length £ asso- 
ciated with an ordered Pfaffian expansion, exactly {2£)\\ are adjacent. 

Proof. To count the total number Na(£) of all adjacent loop- like strings of the length 
£, we consider a specific pair of adjacent loop-like strings of the length £ represented 
by the sequences of arguments 

(zj 1 , Zj 2 ){zj 2 , zj 3 ){zj 3 , Zj 4 ) ■ ■ ■ (zj p l , zj p ) ■ ■ ■ (zj e _ 1 , Zj e )( z ji^ ' %i) 

pair \ 1 pair \\ 2 pair \\ (t—1) 

and 

, Zj 2 ){zj 2 , Zj 3 ){zj 3 , Zj 4 ) ■ ■ ■ (Zj pl , Zj p ) ■ ■ ■ (Zj { _ 1 , Zj e )( z je^ ' z ji)- 
pair i\ 1 pair \\ 2 pair \\ (^-1) 

Here, the mutually distinct j p take the values from 1 to I. The two strings are identical 
up to an exchange of the first and the last arguments Zj t Zj ± . The remaining 2(^—1) 
arguments are distributed between the kernels in such a way that an adjacent string 
is formed in accordance with the Definition 5.7; the (£ — 1) undcrbraces identifying 



f For an example, please refer to the second and third column in (5.11). Also, see Lemma 5.3 for 
the explanation of the number 16 = 4!! 2!! and Definition 5.6 for the notion of a loop-like string. 
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(I — 1) pairs of complex conjugate arguments highlight the structure of an adjacent 
string. 

The total number Na{£) of all adjacent loop- like strings of the length £ equals 
the number of ways to generate those strings from the two depicted above. As soon as 
there are (i) £ ways to assign a number from 1 to £ to the label ji, (ii) {£ — 1)! ways to 
assign the remaining (£ — 1) numbers to the (£ — 1) pairs left (labelled by j 2 , • • • 
and (iii) 2 £_1 ways to exchange the arguments Zj k z± Zj k (k = 2, • • • ,£) within those 
(£ — 1) pairs, we derive: 

Af A (£) =2x/x(M)!x 2 1 - 1 = {2£)\\ (5.24) 
End of proof.J ' 

Remark 5.1. Since the above reasoning holds for loop-like (sub)strings of any length 
1 < p < £, one concludes that the total number of adjacent (sub)strings of the length 
p equals Ma{p) = (2p)H. 

Lemma 5.5. Any loop-like (sub)string of the length p can be transformed into an 
adjacent (sub)string of the same length by means of proper (i) permutation of kernels 
and/or (ii) permutation of intra-kernel arguments. 

Proof. To be coherent with the notations used in the proof of the Lemma 5.3, we 
deal below with a loop-like string of the length £. However, the very same argument 
applies to any loop-like (sub)string of the length 1 < p < t so that our proof (based 
on mathematical induction) holds generally. 

(i) For 1=1 and £ = 2, the Lemma's statement is obviously true. Indeed, a 
loop-like string of the length £ = 1 is automatically an adjacent one. For £ = 2, a 
loop-like string composed of two kernels is reduced to an adjacent string by utmost 
one intra-kernel permutation of arguments. 

(ii) Now we assume the Lemma to hold for loop-like (sub)strings of the length £ 
(that is, that any loop-like string of the length £ can be reduced to an adjacent string 
by means of the two types of allowed operations). 

(iii) Given the previous assumption, we have to show that a loop-like (sub)string 
of the length (£ + 1) can also be reduced to an adjacent (sub)string. It follows from 
the proof of the Lemma 5.3 (see the discussion around (5.22)) that a loop- like string 
of the length (£ + 1) can be generated from a loop-like string of a smaller length £ 
by adding an additional pair of arguments (2^+1, Z£+i) followed by exchange of either 
zi + i or zg + i with one of the arguments of the original loop-like string of the length £. 
Since, under the induction assumption (ii), the latter can be made adjacent, 

(z l+1 ,z e+1 ) <g> [{z h ,...) ■■■ (...,z jq ){z jq ,...) ■■■ (...,z h )}, (5.25) 

S v ' 

adjacent string of the length £ 

one readily concludes that we are only two steps away from forming an adjacent string 
of the length (^+1) out of (5.25). Indeed, an exchange of arguments z e+ i ^ zj q brings 
(5.25) to the form 

{zt+i,z jq ){z h ,...) ••• (...,z t+1 )(z 3q ,...) ••• (...,%) (5.26) 

I In particular, there should exist eight adjacent strings in the ordered Pfaftian expansion for 1 = 2. 
This is in concert with the explicit counting in Example 5.7. 
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which boils down to the required adjacent string 

(z h ,...) ••• {...,z e +i)(z e+1 ,z jq )(z jq ,...) ••• (...,z h ) (5.27) 

V v ' 

adjacent string of the length (^+1) 

upon moving the pair (zt+i, Zj q ) through (q — 1) pairs on the right. Exchanging 
zg + i <^ zj q instead can be done in the same way. ■ 

Remark 5.2. The reduction of a loop-like (sub)string to an adjacent (sub)string 
is not unique. For instance, the loop-like string 

+(13)(23)(12) 

(see the first term in (4.16)) can be reduced, by permutation of kernels and 
permutation of intra-kernel arguments, to one of the following adjacent strings: 

+(13)(32)(21), -(12)(23)(31), 
+(21)(13)(32), -(23)(31)(12), 

+ (32)(21)(13), -(31)(12)(23). (5.28) 

To handle the problem of the non- unique reduction of a loop- like string to an adjacent 
string, the notion of string handedness has to be introduced. 



5.2.5. Handedness of an adjacent substring 

Definition 5.8. Close an adjacent substring S into a loop by "gluing" the right 
argument of the last kernel with the left argument of the first kernel below the chain 
as in Figure 4-., 

S : ~» • (w ni w j2 ) ■ (w h , % 3 ) • • • ,%,)■■■ (Wj^ , w jp ) • (w jp w h ) » ~* 

first kernel last kernel 

(here, a set of the arguments (wi, • • • , W2 P ) is specified by (5.9) with I set to p, and the 
symbol • denotes a gluing point). Read the arguments of a loop, one after the other, in 
a clockwise direction ( as depicted by the symbol ), starting with Wj 1 until you arrive 
at Wj p . If ccr is the number of times an argument Zj q is followed by its conjugate Zj q 
for all q e (1,- •• ,p), 

(• ■ • i Z j q ) ' {Zjq :■■•): 

an adjacent substring S is said to have the handedness B^o^Or), where «l = 
p - a R . 

Example 5.8. The handedness of eight adjacent strings considered in the Exam- 
ple 5.7 with p = I = 2 is listed below: 

H(2,0), H(0,2), H(0,2), H(2,0), 

H(l,l), H(l,l), H(l,l), H(l,l). (5.29) 

Lemma 5.6. Out of (2£)\\ adjacent strings of the length I arising in an ordered 
Pfaffian expansion, there are exactly 

Af a (£) = £1 ( l ) (5.30) 
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strings with the handedness M(a,£ — a). 

Proof. A string with the handedness M(a,£ — a) closed into a loop (see Defini- 
tion 5.8) contains a "left" fragments (. . . , Zj ) • (zj , . . .) and {£ — a) "right" fragments 
(. ..,zj) ■ {zj q ,...) with the opposite order of complex conjugation in the nearest 
neighbouring kernels; each fragment is labelled by an integer number j q G (1, • • • ,£). 
To count a total number of all strings with the handedness W(a,£ — a), we notice that 
there exist ways to distribute "left" and "right" fragments on the loop, and £\ 
ways to assign £ integer numbers (from 1 to £) to the labels j\, ■ ■ ■ Applying the 
combinatorial multiplication rule completes the proof. ■ 

Remark 5.3. It is instructive to realise that the Lemma 5.4 can be seen as a corollary 
to the Lemma 5.6. Indeed, the total number of all adjacent strings of the length £ is 
nothing but 



Not unexpectedly, this result is in concert with the Lemma 5.4. 

Corollary 5.2. The number of adjacent substrings of the length p, 1 < p < £, 
with the handedness H(a,p— a) equals 



The total number of all adjacent substrings of the length p is (2p)\\ 

Proof. Follow the proof of the Lemma 5.6 and the Remark 5.3 with £ replaced 
by p. M 

Example 5.9. Out of 16 longest loop-like strings arising in the Pfaffian expansion for 
£ = 2, there are eight adjacent as explicitly specified in Example 5.7. The handedness 
of those strings was considered in Example 5.8. In accordance with the Lemma 5.6, 
there must exist Ao(2) = 2 strings of the handedness H(0, 2), A/2 (2) = 2 strings of the 
handedness H(2, 0), and A/1 (2) = 4 strings of the handedness H(l, 1). Direct counting 
(5.29) confirms that this is indeed the case. 

5.2.6. Equivalence classes of adjacent (sub)strings 

Having defined a notion of the handedness of an adjacent string, we are back to 
the issue of a non-uniqueness of reduction of a loop-like string to an adjacent one. To 
deal with the indicated non-uniqueness problem, we would like to define, and explic- 
itly identify, all distinct equivalence classes for (2^)!! adjacent strings arising in the 
context of an ordered Pfaffian expansion. 

Definition 5.9. Two adjacent strings, Si and Sj, are said to be equivalent ad- 
jacent strings, Si <~ Sj, if they can be obtained from each other by (i) permutation 
of kernels and/or (ii) intra-kernel permutation of arguments. 





(5.31) 



Example 5.10. The six adjacent strings (5.28) are equivalent to each other. 
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Definition 5.10. A group of equivalent adjacent strings is called the equivalence 
class of adjacent strings. The j-th equivalence class to be denoted ACj consists of 
||-4Cj|| equivalent adjacent strings. 

Lemma 5.7. All (2£)\\ adjacent strings arising in the ordered Pfaffian expansion 
can be assigned to (21— 2)!! equivalence classes {ACi, ■ ■ ■ , AC( 2 £-2)n} °f adjacent 
strings, each class containing 21 equivalent adjacent strings: ||-4Cj|| = 21 for all 

iei,-,(2^-2)!L 

Proof. Let us concentrate on a given adjacent string S with the handedness H(a, £— ex) 
that belongs to an equivalence class ACj and count a number of ways to generate equiv- 
alent strings out of it by means of (i) permutation of kernels and/or (ii) intra- kernel 
permutation of arguments without destroying the adjacency property. Two comple- 
mentary generating mechanisms exist. 

• First way (Ml): One starts with an adjacent string of the handedness H(a,t — a), 

~* • ( w h > w h ) • (w j2 ,%,)■■■ (wj^ ,wj p )--- (w Je _ 17 w H ) ■ (w jt , w h ) • ~» (5.32) 

S v ' S v ' 

first kernel last kernel 

(the reader is referred to the Definition 5.8 for the notation used), simultaneously 
flips the intra-kernel arguments in all £ kernels, 

{W h , Wjl ) ■ (% 3 , W j2 ) ■ ■ ■ (w jp , W Jp _ t ) ■ ■ ■ (W je , Wj^) ■ {w h , w jt ) , 
first kernel last kernel 

and further permutes the kernels in a fan-like way so that the last ^-th kernel 
in the string becomes the first, the (£ — l)-th kernel in the string becomes the 
second, etc.: 

~* • iw h , w jt ) ■ (w je , Wj^J ■ ■ ■ (w 3p , w jp _ 1 ) ■ ■ ■ (w j3 , w h ) • (w n , w n ) • ~* (5.33) 

S v ' S v ' 

last kernel first kernel 

The so obtained adjacent string is equivalent to the initial one (5.32) but possesses 
the complementary handedness M(£ — a, a). 



• Second way (M2): One starts with an adjacent string of the handedness H(o;, £ — 
a), and permutes the £ kernels in a cyclic manner to generate (£ — 1) additional 
(but equivalent) adjacent strings with the same handedness (to visualise the 
process, one may think of moving a gluing point • through the kernels in (5.32)). 

The two mechanisms, Ml and M2, combined together bring up 2£ equivalent adja- 
cent strings due to the combinatorial multiplication rule. As a result, wc conclude 
that H-ACj II = 21. Consequently the number of distinct equivalence classes equals 
{2l)\\/(2£) = (2l-2)\\. ■ 

Example 5.11. To illustrate the Lemma 5.7, we consider an ordered Pfaffian ex- 
pansion for £ — 2 as detailed in the Example 5.1. The eight adjacent strings of the 
length £ — 2 were specified in the Example 5.7; their handedness was considered in 
the Example 5.8. In accordance with the Lemma 5.7, there should exist 2 distinct 
equivalence classes, each containing 4 adjacent strings. Indeed, one readily verifies 
that those two equivalence classes are 

Ad: -(12)(21), -(21)(12), -(T2)(21), -(21)(12) 
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Longest Loop-Like Strings 

Total amount: (2£)\l(2£ - 2)!! 



Adjacent Strings 

Total amount: (2£)H 



u 



Non- Adjacent 
Strings 



r ci 



l «-(2<?-2)!! ) 



Ad 



AC 



£ strings: H(ai,£ — a?i) 
£ strings: H(£ — ai,ai) 



£ strings: M(aj,£ — ctj) 
£ strings: M(£ — otj,otj) 



(2£-2)\ 



£ strings: M(a (2 i-2y.'.J - a^e-2y.\) 
I strings: U(£ - a^-2)n, a(2*-2)n) 



Figure 5. Any longest loop-like string from the equivalence classes 
{Ci, • • • ,C(2i— 2)!l} can t" 3 reduced to an adjacent string belonging to one 
of the equivalence classes {AC\ , ■ ■ ■ , AC^e-2)'.'. } °f adjacent strings by means 
of proper permutation of kernels and/or permutation of intra- kernel arguments 
(Lemma 5.5). This reduction is non-unique (Remark 5.2). Each equivalence class 
ACj of adjacent strings contains i adjacent strings of the handedness W(otj,£ — <Xj) 
and £ adjacent strings of the complementary handedness H(£ — ctj, ctj) (Remark 
5.4). Notice that two distinct equivalence classes AC{ and ACj (where i ^ j) 
may have the same values for a: on = ctj . 



and 

AC 2 : +(12)(21), +(21)(12), +(12)(21), +(21)(12). 

Remark 5.4. In fact, a generic prescription can be given to build (21 — 2)!! dis- 
tinct equivalence classes {AC\, ■ ■ ■ , AC( 2 i-2)\\} f° r (2^)!! adjacent strings arising in the 
ordered Pfaffian expansion. Because of the "duality" between equivalent adjacent 
strings with complementary handedness H(a, i — a) and H(£ — a, a) discussed in the 
proof of the Lemma 5.6, the adjacent strings whose handedness M(a, I— a) is restricted 
by the inequality < a < [£/2\ will form a natural basis in the consideration to follow. 



• The case £ = 2A + 1 odd. 

(i) For < ai < A, (i.a) pick up an adjacent string with the handedness 
M(ai,£ — a\), out of A/" Ql (£), and generate I equivalent adjacent strings with 
the same handedness through the mechanism M2 of the Lemma 5.7. (i.b) Apply 
the mechanism Ml of the same Lemma to each of the £ adjacent strings generated 
in (i.a) to create I more equivalent adjacent strings of the handedness W(l—a\,a{) 
hereby raising their total amount to 21. The strings generated in (i.a) and (i.b) 
are said to belong to the equivalence class AC\. 
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(ii) To generate the next equivalence class AC2 ^ ACi, pick up an adjacent string 
not belonging to AC\ with the handedness H(a2, 1 — 012) out of (A/ Q1 (t) — i) left 
(again, a 2 is restricted to < a 2 < A), and repeat the actions described in (i.a) 
and (i.b) to generate another set of 21 equivalent strings. These will belong to 
the equivalence class AC2 which is distinct from AC\. 

(iii) To generate the j-th equivalence class ACj, one picks up an adjacent string 
out of (J\f ai {£) — (j — 1)^) left an d repeats the actions sketched in (ii). 

(iv) The procedure stops once there are no adjacent strings left. Obviously, the 
total number of equivalence classes is 

1 11/21 A /2A + 1\ 

7 = (2A)l£( 2 = (2A)!2 2A = (2^-2)!! (5.34) 

This is in concert with the Lemma 5.7. 



• The case I = 2A even. 



In this case, special care should be exercised for the set of adjacent strings with 
the handedness H(A, A) because these adjacent strings are seZf complementary: 
the mechanism Ml applied to any of those adjacent strings generates a string 
with the same, not complementary, handedness. The latter circumstance can 
readily be accommodated when giving a prescription for building (2£ — 2)!! dis- 
tinct equivalence classes of adjacent strings. 

(i) First, we separate all adjacent strings with the handedness W(a,£ — a) where 
< a < A — 1 and apply a procedure identically equivalent to that described for 
the case I odd to generate distinct equivalence classes of adjacent strings. The 
total amount of distinct equivalence classes built in this way equals 

U = ) 1 E _ \ a (t) = (2A - 1)! g = (2t 2)!! - . (5.35) 

q=0 a=0 ^ ' ^ ' 



(ii) Second, having generated in the previous step L\ distinct equivalence classes 
{ylCi, • • • , ACl-l} of adjacent strings, we concentrate on the adjacent strings with 
the handedness H(A, A) not treated so far. To this end, we (ii.a) pick up an adja- 
cent string out of A/a (2 A) with the handedness H(A, A) and perform the operations 
Ml and M2 to generate 21 = 4A equivalent strings with the same handedness. 
The 21 equivalent adjacent strings will belong to a certain equivalence class, say, 
ACl 1 +i- (ii-b) In the next step, we pick up an adjacent string with the handed- 
ness H(A, A) out of (A/a(2A) — 21) left, and perform the operations detailed in (ii.a) 
in order to generate yet another set of 21 equivalent adjacent strings belonging to 
an equivalence class ACl x +2- (ii-c) We proceed further on until the last available 
equivalence class composed of 21 adjacent strings is formed, ^4Cl 1+ l 2 , where L 2 
equals 

U _l mti .Q^pX\. (5.36) 



21 AW 2 V A 
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Hence, for £ = 2A, the total number of equivalence classes for adjacent strings 
equals 

L 1 + L 2 = (2£-2)\\ (5.37) 
as expected from the Lemma 5.7. 

5.3. Integrating out all longest loop-like strings of the length £ 

More spadework is needed to prove the Pfaffian integration theorem. Below, we will 
be interested in calculating the contribution Cl(£) of longest loop-like strings (of the 
length £) into the sought integral (5.2b): 




Qn(Zi, ZjJ Qn(zi 7 Zj) 
Qn (%i : Zj ) Qn(z% > Zj ) 



11-kUJ longest 

loop — like strings 



-1 



]T sgn(a) [] Q 

n (w a ( 2 j+l), «V( 2 j+2)). (5.38) 

ves' 2e j=o 

In the second line of (5.38), only that part of the ordered Pfaffian expansion (5.8) 
appears which corresponds to a set of all loop-like strings of the length £. They are 
accounted for by picking up proper permutations S' 2e C S^t in the expansion (5.8), 

S'21 1 — * (longest loop — like strings of the length £) . 

Although, in accordance with the Lemma 5.3, the number of terms in the expansion 
(5.38) equals j\f L (£) = (2£)\\(2£ - 2)!!, there is no need to integrate all of them out 
because various loop-like strings belonging to the same equivalence class yield identical 
contributions. The latter observation effectively reduces the number of terms in (5.38) 
so that 

(2£)" f e e i 

Zt - JC 3 = 1 aeS' 2 ' e 3=0 

Here, the prefactor (2£)\\ = 2 l £\ equals the number of longest loop-like strings in each 
equivalence class; the cr-series runs over the permutations S 2l C S' 2l corresponding to 
Ul{£) = (2£ — 2)!! longest loop-like strings, each of them being a representative of one 
distinct equivalence class, 

&2t ^ € Ci, ■ ■ ■ ,<S(2£_2)M € C( 2 ^_2)!l}, (5.40) 

see the Lemma 5.3. There are (2£ — 2)!! terms in (5.39). 

To perform the integration explicitly, one has to reduce the longest loop- 
like strings in (5.39) to the form of adjacent strings as discussed in the Lemma 
5.5. In accordance with the Lemma 5.7, there exist (21 — 2)!! equivalence classes 
{ACi, ■ ■ ■ ,ACr2i-2)u} OI adjacent strings, each of them containing 2£ equivalent 
adjacent strings (see also Fig. 5). This results in the representation 
„ t i-i 

C L (£)= I JJ dTT(^) Yl S S n ( £T ) YlQn(Wa(2j + l),W^ 2j+ 2)), (5.41) 
JC 3 = 1 aeS- e 3=0 

where the cr-series runs over the permutations S' 2e C S' 2t corresponding to {21 — 2)!! 
adjacent loop-like strings of the length £, each of them being a representative of each 
one of existing equivalence classes of adjacent strings, 

S' 2 ' t i ^ {Si e AC u ---,S {2 i-2)u G AC (2e - m }- (5.42) 
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The number of terms in (5.41) is (21 - 2)!!. 

To proceed, we rewrite (5.41) in a more symmetric form that treats all adjacent 
strings on the same footing: 

1 r 1 - e ~ 1 

CL ^ = 2l / 1^ d7T ( Z ^ E Sgn ( CT ) Y[Qn(Wa(2j+l),W rT{ 2 ] +2))- (5-43) 

JC J=l <t£S'& j=0 

Here, the cr-series runs over the permutations S' 2l C S 2e corresponding to all adjacent 
loop-like strings of the length t. 



1 2( 



{{S?\ • • • , sg>} e AC,, • • • , { sr-V»\ ■ ■ -A ( r m) } e AC {2e . 2) ,} . (5.44) 

As soon as there exist 21 equivalent adjacent strings in each equivalent class ACi 
of adjacent strings, the prefactor {2i)~ 1 was included into (5.43) to avoid the 
overcounting. 

An advantage of the representation (5.43) can be appreciated with the help of the 
Lemma 5.6. According to it, the summation over the permutations a G S 2 " can be 
replaced with the summation over all longest adjacent strings with a given handedness 
M(a,e-a), for all a E (0,f): 

E S § n ( CT ) HQn(Wa(2 ] + l),W a{2:j+2) ) = Y E ( 5 ' 45 ) 
ceS' 2 'l 3=0 a=0 i=l 

Here, <Sj(a) denotes the i-th adjacent string of the length ||<S,(a)|| = I with the 
handedness M(a,£ — a). In accordance with the Lemma 5.6, 

jV a (e) = e\( e ). (5.46) 



Given (5.45), the integration in (5.43) can be performed explicitly. Due to the new 
representation 

i jV a (i) e 

ClW = Yf E E / II Mzi)Si(a), (5.47) 

a =0 i=l Jc j = l 

one has to calculate the contribution of a string Si(a) with the handedness H(a, t — a) 
to the integral: 

t 

h(a) = / TT dk{zj)Si{a). (5.48) 

(i) The case a — is the simplest one. Having in mind the definition (5.2a) and 
introducing an auxiliary matrix <f with the entries 

9* = o / dlK ( z )0j( z )Qk(z), (5.49) 
L Jc 

we straightforwardly derive: 

~e-i 



1/(0) = sgn(<7 ) / f[dir{zj) 

JC 3 = 1 



Y[ Qn{zj,z J+1 ) Q n {z e , zi) = -tr ( o )Tl _i) [(A^) £ ] -(5.50) 

_3 = 1 

Here, the permutation sign, sign(cro), is sgn(ero) = — 1 (see (5.9)), while the trace 
tr(o, n -i)(' • •) reflects the fact that the integrated adjacent string is loop-like. Impor- 
tantly, the result of the integration (5.50) does not depend on a particular arrangement 
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of the arguments zj and Zj as far as the handedness H(0,£) is kept. 



(ii) The case a = i associated with the adjacent loopdike strings of the handedness 
M(£, 0) can be treated along the same lines to bring 



(5.51) 



(iii) More care should be exercised for < a < £. In this case, two adjacent strings with 
the same handedness H(a, I — a) may bring different contributions into the integral 
(5.48). For instance, the adjacent string § 



Si (a) = sgn(a a ) 



j=a+l 



Y\ Qn(Zj, z j+l) Qn{z a ,Z a + i) 
J = l 

yields the contribution 

if\ a ) = (-lr-hr^^ [(£<r)<w- a ] . 

At the same time, the adjacent string 
S ( e a )( a ) = sgn(<7 a ) 



Qn{z£, Zi) 



(5.52) 



-a-1 



II Qn(Zj,Zj+l) 
3 = 1 



Qn{zt-a, Z£- a+ i) 



1 i Qn{Zj 
j=£-a+l 



Qn{zt, Z\), 



possessing the very same handedness, yields 

IP\a) = (-ir-h Ho ,n-i) [(£?) /_q (£O q ] • (5.53) 

In deriving the results (5.52) and (5.53), we have used the fact that the permutation 
sign sign(cr Q ) is sign(er Q ) = (-1) Q_1 . 

It can readily be seen that the adjacent strings with a given handedness H(a, £— a) 
bring all possible (^) contributions, or words, that can be represented as a trace 

Wj(£,a) = -tr (0>n _i) [• • • (-AT) • • • (AO • • • (-AO • • •] (5-54) 

N ^ ' 

a letters ( — p,s*) and £ — a letters (/}<?) 

of a product of a matrices, or letters, (— A<?*) and £ — a matrices (letters) (p,<H) 
distributed in all possible (*) ways. Hence, the index j in (5.54) takes the values 
1 < j < (^). Importantly, each word Wj(£,a) appears exactly i\ times since there 
always exist l\ adjacent strings Si k (a) (k = 1, • • • , £\) which are related to each other by 
a permutation of the integration variables in (5.47). As a result, the latter is reduced 
to 



(i) 



a=0 j=l 



Spotting that 



a) 



W j( £ > a ) = - tr (0,n-l) {flS-flC) 



(5.55) 



(5.56) 



The empty products are interpreted to be 1. 
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we obtain 

C L (£) = -\(t- 1)! <*((>,„_!) [(A? - AT)'] ■ (5.57) 

Finally, noticing that the matrix — jii<f*) under the sign of trace is related to the 
matrix v defined by (5.2c) as 

t> = 2i(£?-AO. ( 5 - 58 ) 
we arrive at the remarkably compact result 

C L (£) = -l^0£tr {o , n _ 1) (v e ). (5.59) 
Hence, we have proven the following Theorem. 

Theorem 5.2. Let dn(z) be any benign measure on z e C, and the function Q n (x, y) 
be an antisymmetric function of the form 

^ n— 1 

Qn{x,y) = 2 XI Qj( x )^jkqk{y) 
j,k=o 

where qj 's are arbitrary polynomials of j-th order, and p, is an antisymmetric matrix. 
Then the integration formula 



c L (£) = n / d < z ^ ( p f 

j=1 Jzjec \ 



l 



2Cx2t J longest 

loop — like strings 



tT^.n-l)^) (5-60) 



2 (2z) 

holds, provided the integrals on the l.h.s. exist. Here, the matrix v is determined by 
the entries 

n-l 
k=0 



n—1 

i^Aa.fc / dir(z) [q k (z)qp{z) - qp(z) q k (z)\ . 



The following corollary holds. 

Corollary 5.3. Consider a set of 1p arguments 

{wi = zi, w 2 = zi, ■ ■ ■ , w 2p -i = z p , w 2p = z p }, (5.61) 

all (2p)\ permutations of which are denoted by S 2p . Take the subset S' 2p C S 2p of S 2p 

corresponding to all loop-like strings of the length \\S^\\ =p, 
p-i 

S i P) = J\Qn{w a . {2 j +1) ,W a . (2j+2) ). 

Here, Oi labels the i-th permutation cr, € S' 2p - The Theorem 5.2 implies: 
I r p P- 1 

Cl< &) = ^p—{ I II dll ^ Z ^ H Sgn ( CT ) II Qn{Wc( 2 j+l),W a(2i+2) ) 

P ' Jc i=i °es> 2p j=o 

l(v- 1)' 

= -2^-tr (0;n _ 1) (^). (5.62) 
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5.4- Integrating out compound strings 

Having dealt with the longest loop-like strings arising in the ordered Pfaffian expansion 
(see Theorem 5.2) and having extended the result of this theorem to loop-like strings 
of a smaller length (see Corollary 5.3), we now turn to the treatment of the remaining 
(2£)\ — {2£)\\{2£ — 2)!! strings to be referred to as compound strings. Our final goal 
here is to calculate their contribution to the integral in the l.h.s. of (5.26). 

To give a definition of a compound string, wc remind that any given string 5, of 
the length = I from the ordered Pfaffian expansion can be decomposed into a set 
of loop-like substrings of smaller lengths (see Lemma 5.2). According to the Lemma 
5.5, each of the above loop-like substrings <S> 4 can further be reduced to the form of 
an adjacent substring. This leads us to the following definition: 

Definition 5.11. A string Si of the length \\Si\\ — £ from the ordered Pfaffian expan- 
sion is called a compound string if it is composed of a set of adjacent substrings sf ] ^ 
of respective lengths \\sf 3 \\ = £j such that Si = U . sf 3 ^ with J2j &j = 

Remark 5.5. (i) This definition suggests that all compound strings can be clas- 
sified in accordance with all possible patterns of unordered partitions A of the size 
|A| = £ of an integer £: 

A=(^V,^ 9 ). (5.63) 

The frequency representation (5.63) of the partition A says that the part £j appears 
Oj times so that 

£ = f^£ j a j . (5.64) 
j=i 

Here, g is the number of inequivalent parts of the partition A. 

(ii) Alternatively, the partition (5.63) of an integer £ can be represented as 

A = {£ u ■■■Jr), (5.65) 

where the order of Ij's is irrelevant, and some of them can be equal to each other. 
Obviously, 

r 
3=1 

Here, r is the length of unordered partition A. 

(iii) The correspondence between compound strings of the length £ and unordered 
partitions A of the size |A| = £ gives rise to a topological interpretation of compound 
strings, which can conveniently be represented in a diagrammatic form (see Fig. 6). 
The diagram for a generic compound string 

a IK- X> '• 

3=1 3=1 

consists of r loops, the j-th loop depicting an adjacent substring Sj of the length £j. 
Such a diagram will be said to belong to a topology class {£\, ■ ■ ■ ,£ r } = , • • • , £g 3 }- 
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0.0 




Figure 6. A diagrammatic representation of a compound string belonging to the 
topology class A = (f^ 1 , • • • ,1-g 3 ) with X)f=i = ^ an< ^ 5Zf=i °j = r - Here, 
<; denotes the number of inequivalent parts of the partition A whilst r equals the 
total number of loops. 



Remark 5.6. The number of tope-logically distinct diagrams equals the number 
p{£) of unordered partitions of an integer t (see Remark 5.5). These are known to 
follow the sequence 

p{£) = {1, 2, 3, 5, 7, 11, 15, 22, 30, 42, • • •}. (5.67) 

An exact evaluation of p{£) can be performed with the help of Eulcr's generating 
function (Andrews 1998) 

oo oo 1 

£^=1171— ( 5 - 68 ) 

1=0 m=l [l q > W°° 

The asymptotic behaviour of p{£) for I » 1 was studied by Hardy and Ramanujan 
(1918),' 

m ~ ST! cxp ( 7rv/ ^) • (5 - 69) 

To calculate the contribution of compound strings to the integral in the l.h.s. of 
(5.26), one has to determine (i) the contribution of a diagram belonging to a given 
topology class to the integral, and (ii) the number of diagrams within a given topology 
class. 

Lemma 5.8. The number of diagrams belonging to the topology class X = 
(£l\---,£g 9 ) equals 

Proof. To determine the number Af^ a i ...fn} of diagrams belonging to a given 
topology class ■ ■ ■ ,£ r } = {£°^ , ■ ■ ■ 7 £g°}, we use the multiplication principle. 

• First, we distribute the pairs of arguments {£~i,- ■ ■ ,£ r } between r loops. This can 
be achieved by mi ways, 

m, = ( I )('- h ) ir - 2 ) (<-*-"- H x I. 



2 / V «r-l 



£ r J r\ 
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The factor 1/rl reflects the fact that the order of the r loops is irrelevant. (Indeed, 
a given topology class is associated with an unordered partition of an integer £) . 
Simple rearrangements show that m\ simplifies down to 

mi = i^ri7v (5J1) 

• Next, we shuffle all inequivalent loops (those with distinct lengths). This can be 
achieved by m 2 ways, 

(a 1 + --- + a g )\ r! 
m 2 = : : = — : r. (5.72) 

(J\\---0-g\ a\\---CTg\ 

As a result, the total number of diagrams belonging to the topology class A = 
equals 

= ™ = iii e [ llail 1 .. (Tg r (5-73) 
The observation ]J r IA = JlLi^'D^ cnds thc P roof - ■ 



Lemma 5.9. A diagram associated with a topology class A = (l^ 1 , ■ ■ ■ , £g e ) con- 
tributes 

to the integral in the l.h.s. of (5.2b). The function Cl{p) is defined by (5.62) of the 
Corollary 5.3. 

Proof. The above claim is a direct consequence of the Lemma 5.2, Definition 5.11, 
and Remark 5.5. End of proof. ■ 



Theorem 5.3. All compound strings belonging to the topology class A = {il 1 ,■■■,£', 
yield, after the integration in (5.2b), the contribution 



1 / 1 r~l-s 

-7T tr (o,«-i)(^ J ) 



f/aA\ 2 



(5.75) 



Proof. Observe that C\ = ... fgyC^i ... ^sj, and make use of (5.62), (5.70) 

and (5.74) to derive (5.75). End of proof. ■ 



5. 5. Proof of the Pfaffian integration theorem 

Now we have all ingredients needed to complete the proof of the Pfaffian integration 
theorem announced in Section 5.1. Indeed, in accordance with the topological 
interpretation of the terms arising in the Pfaffian expansion (5.8), the integral on 
the l.h.s. of (5.26) is given by the sum of contributions of adjacent strings of two 
types: the longest adjacent strings and the compound strings. The contribution of 
the former, Cl{£), is given by thc Theorem 5.2 while thc contribution of the latter, 
C\, is determined by the Theorem 5.3. As soon as 

Cl{£) = Ca=(*i), 
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one immediately concludes that the integrated PfafBan equals the sum of C\ over all 
unordered partitions A of the size | A| = I: 



n rf7r<z j) pf 



i=i 



Qn {%i 7 Zj ) Qn > Zj ) 
Qn (-^i j Zj ) 7 Zj ) 



|A|=£ 



e n 

|A|=0'=1 



2^ 



tr 



(0,n-l) 



(<*) 



(5.76) 



Quite remarkably, the r.h.s. of (5.76) can be recognised to be a zonal polynomial 
(Macdonald 1998) 

\\\=e j=i 

with the arguments 



1 



Pi = 2 tr (o,n-i)(^) ; i = l,2,- 



As a result, we conclude that 



,1 1 i 

■^(i") ( o tr (o,n-i)W , • • • , -tr ( o, n _i)t; 



This coincides with the statement (5.26) of the Pfaffian integration theorem. 



(5.77) 



(5.78) 



6. Probability function p n ,k'~ General solution and generating function 

6.1. General solution 

The general solution for the probability function p n ^ of a fluctuating number of real 
eigenvalues in spectra of GinOE is now straightforward to derive. Indeed, it was shown 
in Section 3.3 that the probability function p Ut k admits the representation 



Pn,k 



n 



d Zj 



j—l Jim. Zj >0 



x erfc 



iV2 



pf 



T^n {%i i Zj ) T-^n i Zj ) 
T^n i^Zi , Zj ) T^n {. Zi , Zj ) 



(6.1) 



2ix2£ 



with the kernel function T> n (x,y) given by (3.12) and (3.13). 

The ^-fold integral in (6.1) can explicitly be performed by virtue of the Pfaffian 
integration theorem after the identification 

dir(z) = e- (z2+22)/2 crfc e(lmz)d 2 z, Q n {x,y) = e^ 2+y2 ^ 2 V n {x,y). 



V «V2 

Straightforward calculations bring 



Pn,n r, 
Pn,k = A 



(1<) 



-tr 



(O.n-l) 1 



-tr 



(O.ra-l) 1 



(6.2) 



where the matrix v is given by (5.2c). Combining the definition (5.2c) with (3.21), 
(3.22), (3.23), (3.25) and (3.27), one concludes that 



v = 2i{jxx) = a-, 



(6.3) 
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see Appendices A and B. Finally, making use of the trace identity 

tr^-i)^ = 2tr (0;L „ /2J _ 1) £ j (6.4) 
proven in Appendix C, we end up with the exact formula 

Vn.k = ^tr ( 0) [n/2j -l) Q 1 ! ' ' ' i tr (0, L"/2J -1) ^) ■ (6-5) 

The entries of the matrix q, calculated in Appendix C, are given by 

Q e JT= / dyy 2 ^- 1 ^^ eric (yV2) 
Jo 

x "(2a + 1) L^ + T ) - 1 (-2 2 / 2 ) + 2y 2 L^_T )+1 (-2» 2 )] (6.6) 

and 

-odd _ -even ( A\m-0 TO ' -oven / fi y\ 

Q a ,p-Q a ,f3 \ V (2m)! Q ' m ^ 

for n = 2m even and n = 2m + 1 odd, respectively. 

t>. 2. Generating function for p n ^ 

Interestingly, the entire generating function 

L«/2J 

G n (z) = ze Pn,n-2l (6.8) 

e=o 

for the probabilities p n ^ can explicitly be determined. To proceed, we make use of 
the summation formula 



oo 

Z 



E h ■■■,Pr) = exp I £(-1)^^- I (6-9) 



r 

r = 



well known in the theory of symmetric functions (Macdonald 1998). With p r 
tr (o,|_n/2j-i)f? r j the r.h.s. transforms to 



CXP tr (0iL „/ 2 j-l) E^ 1 )^ 1 ^ - = ° XP ( tr (0,Ln/2j-l) log(l|n/2J + 
' r>l ' 



= det [i + ^] L „ /2JxLn/2J , (6.10) 
resulting in an amazingly simple answer: 

in/2] 

Gn(z) = Z Wn-2e = Pn,n det [1 + *$] [rl/2j x [n/2j . (6.11) 



6. 3. Integer moments of the number of real eigenvalues 

The result (6.11) allows us to formally determine any integer moment E[A/"*] of the 
fluctuating number j\f r of real eigenvalues in the spectra of GinOE. Denoting the 
fluctuating number of complex eigenvalues through 2jV c , we derive: 

E[AA«] = E[(n - 2jV c ) q ] = ^ ( q \ q - j {-2)m[M 3 c \. (6.12) 

3=0 
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Since 



dz J 

the formula (6.12) simplifies to 

nK q ] =Pn,n (n-2z^j det [i + ^] Ln/2JxL „ /2J \ z=i - (6.14) 
Of course, to make the formulae (6.11) and (6.14) explicit, the determinant 

d n (z) = det [i + *e] Ln/2JxLn/2J 

has to be evaluated in a closed form. There are a few indications that this is a 
formidable task, but we have not succeeded in the calculation of d n (z) yet. 

7. Asymptotic analysis of the probability p n >n _ 2 to find exactly one pair of 
complex conjugated eigenvalues 

To determine a qualitative behaviour of the probability function p n ,k, an asymptotic 
analysis of the exact solution (1.1a) is needed. In this section, the simplest probability 
function p„,„- 2 , 

/>oo 

Pn.n-2 = 2p n ,n / dy y e y ' erfc(y72) L 2 n _ 2 (-2y 2 ), (7.1) 
Jo 

is studied in the largc-n limit. Our consideration is based on an alternative exact 
representation for p n , n -2 (see the Theorem 7.1) which is more suitable for obtaining 
regular large-n asymptotics. 

7.1. Alternative exact representation o/p„ ; „_ 2 
Let us define the sequence 

POO 

S n = / dyy e y 2 erfc(yV2)L 2 n (-2y 2 ), n = 0,l,---, (7.2) 
Jo 

such that 

Pn,n-2 = 2Pn,n S n - 2 . (7.3) 

To find an exact alternative representation for S n (and, hence, for p n ,n-2), we (i) 
introduce a generating function t(z) in the form 

oo 
n=0 

which is supposed to exist in some domain Q T e K of the real line R (to be specified 
later on), (ii) calculate t(z) explicitly, and (iii) expand it back in z e Q T . 

Lemma 7.1. The generating function t(z) reads: 



wftere 

-Kz<i (7.6) 
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Proof. The identity || 

J2 L X n {x) z n = (1 - z)-*- 1 exp , \z\ < 1, 

applied in the context of (7.2) and (7.4), gives rise to the representation 

.l + z 



1 f°° ( 1 

t{z) = JYZ z )3 J dyycric(yV2) exp ^y 2 — 



z| < 1. (7.7) 



A change of the integration variable y to 



2 l + z 



i = y>— (7.8) 



followed by integration by parts results in 

T{Z) = 2(l-z*)(l-z) 
where 



5 erfc (a z VT) - y ^ ^ crfc (a 2 



(7.9) 



1 - z 



l + z 

For the boundary term in (7.9) to nullify at £ = oo, the parameter z has to belong to 
the domain 

n T : -l<z<\. (7.10) 
Performing the remaining integral, we end up with (7.5). End of proof. I 

Having determined the generating function t(z), we are going to Taylor-expand it 
around z = in order to arrive at an alternative formula for the sequence S n . As t(z) 
is a relatively simple function, we may expect that S n obtained in this way will also 
have a relatively simple form. 

Lemma 7.2. The following formula holds: 

1 L " /2J / 2 \ 1 / x 

Sn = ^ £ 3^-/ 2 /W (^) 40/2J+1). (7.11) 

Here, a n = \n/2] - [n/2\ . 

Proof. To expand the function t(z) given by (7.5) around z — 0, we represent it 
in the form 



t(z) = ±t 1 (z)(-1 + V2t 2 (z)), (7.12) 



where 



2^ \ > ^) = 1 7^4, (7-13) 



(l-z 2 )(l-z)' Vl-3z 
and constantly use a variant of the Cauchy formula 



J2 a kz k I • [J2hz k I = £c„z", c„ = £a fe &„- fc , ( 7 - 14 ) 



\fe=0 / \fe=0 / n=0 k=0 

See Eq. (5.11.2.1) in Prudnikov, Brychkov and Marichev (1986). 
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where absolute convergence of the resulting series is assumed. 
Expansion ofr\{z). To determine the coefficients cj^ in the expansion 

oo 
k=0 

we notice that 

oo oo 
k=0 k=0 

so that, in the notation of (7.14), 

l + (-l) fc 
a fe = 1, b k = . 

Straightforward application of the Cauchy formula yields 



3=0 



1 

- j ~ 2 



k + l + 



1 + (-!)* 



(7.15) 



(7.16) 



(7.17) 



Lfc/2J + 1. (7.18) 



Expansion of r^(z). To determine the coefficients in the expansion 

oo 

(2)„fc 



we notice that 



fe=o v 7 
so that, in the notation of (7.14), 



==E(~!/ 2 )(- 3 *) fc 



■1/2 



The Cauchy formula yields 

„(2) _ ' 



c fe = Yl a i bk -i = 3 



fc(V2)fc p / 1 ,1 , 1 



(7.19) 



(7.20) 



(7.21) 



(7.22) 



The latter can be expressed in terms of Legendre polynomials by means of the 
identity % 



2Fl (- l r~ k -\- k > w )- (i/2 )k 

that simplifies (7.22) to 

2 



w 



k/2 



f w + l\ _ fw + l\ 



i" = 3" ! fl|-ij-j"-"' 1 P,- Hii _ 



(7.23) 



(7.24) 



1 Sec Eq. (7.3.1.153) in Prudnikov, Brychkov and Marichcv (1990) 
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Expansion of the product t 1 {z)t 2 {z). To determine the coefficients cj? in the ex- 
pansion 



Tl (z)r 2 (z) = J2c^ z\ 

k=0 

we again use the Cauchy formula 

k 



4 3) = £4 : 

3=0 



(1) J2) 
J J ' 



(7.25) 



(7.26) 



with Cj. and Cj. given by (7.18) and (7.24), respectively. Lengthy but straightforward 
calculations result in 



Lfe/2j 

4 3) = £3^-/ 2 p 2j . +c 

where a fe = [fc/2] - [fc/2j. 



(*) 



(7.27) 



The observation 



S - — c< 3 > - - c (1) 



(7.28) 



completes the proof. 



Theorem 7.1. The probability p n , n -2 to find exactly one pair of complex conjugated 
eigenvalues in spectra of GinOE admits the following exact representation: 



Vn.ri—2 Pn,r< 



L"/2J-1 / „ v 



n/2j 



(7.29) 



Here, a n — \n/2] — [n/2\, and P n stands for Legendre polynomials. 
Proof. Use the Lemma 7.2 and relation (7.3) to deduce (7.29). ■ 
7.2. Asymptotic analysis o/p„ ; „_2 

The result (7.29), combined with the integral representation of Legendre polynomials 

p n {<t>) = ^J de (4> + ^4> 2 - l cose)", (7.30) 

is particularly useful for carrying out an asymptotic analysis of the probability p n , n -2 
in the large-n limit. Indeed, (7.30) facilitates performing a summation in (7.29) leading 
to 

~V2 r +1 



Pn,n — 2 Pn,7i 



dx 



[x + 2) a " -(x + 2) n 
1 - (x + 2) 2 



[n/2\ 



(7.31) 



The large-n behaviour of the integral 

r+1 dx (x + 2) Q " -(x + 2)" 
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is of our primary interest. A saddle-point analysis shows that the saddle point 
x sp ~ —2 lies away from the integration domain x € (— 1,+1). As a result, the 
contribution of the end points of the integration domain, = — 1 and ir = +1, 
should be examined. One can see that the close vicinity x = 1 — e of ir = +1 
dominates exponentially in n > 1. Indeed, the vicinity e G (0, c ) yields 

r° de 3"(l-6/3)"-3""(l-e/3) a " 
""Jo v^ 3 ^) 32(l-c/3)2-l 

where Co is a proper cut-off. In the large-n limit, only a region of order n _1 , e = r/n, 
effectively contributes J„ reducing it to 

3«+l/2 



3" f 00 dr 
\/2n Jo Vt 



r/3 

Combined with (7.31) and (7.32), this estimate leads to the following theorem. 



(7.33) 



Theorem 7.2. The leading large-n behaviour of the probability p n ,n-i is given by 
the formula 

371+1/2 

Pn,n-2 ~ . Pn,n (7.34) 

8 -1/7772 



where p n , n = 2-"(™- 1 )/ 4 . 

Remark 7.1. The Theorem 7.2 implies the inequality p n , n -2 ^> Pn,n- 



8. Correlations of complex eigenvalues of a matrix without real 
eigenvalues 

8.1. GOE correlations in GinOE spectra 

One of the earliest results on eigenlevel statistics in GinOE is due to Ginibrc (1965) 
who spotted that spectra of random real matrices which happened to have no complex 
eigenvalues exhibit the famous GOE behaviour. Indeed, for 7i n e T(n/n), the j.p.d.f. 
(2.5) reduces to + 

0-n(n+l)/4 ™ ™ 

^ n (A 1 ,---,A w )= I] lAi-Ajinexpt-A^). (8.1) 

• 1 Aj = l U/ > i>j = l j = l 

The GOE spectral correlations readily follow (Mehta 2004). 

8.2. GinSE-like correlations in GinOE spectra 

Below, we concentrate on just the opposite case of random real matrices Ho € T(n/0) 
whose spectrum occasionally contains no real eigenvalues. The j.p.d.f. of all complex 

+ As a side remark, we notice that the explicit formula (3.2) for the probability p n ,n can easily be 
derived by integrating P / H n over all of its arguments. Due to Selberg's integral (Mehta 2004) 

f[ f d\ 3 e~ A ?/ 2 f[ | A, - A, | = 2"/ 2 n! f[ T(j/2), 

j = l jR l >j = l j = l 

one obtains p n , n = 2-™("~ 1 )/ 4 . 
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eigenvalues of Hq can also be deduced from (2.5), the result being 

2 -«(n-l)/4 -n/2 J_ 
P-H (Z1, ■ ■ ■ ,Zl) = — — rfn — J } 



(n/2)!m=ir(j/2) 



Zi Zj Zj Zj 



3=1 



Zj — Zj) erfc 



i>j=l 

iV2 



cxp - 



Z 3+*3 



(8.2) 



with n even, n = 2£ Remarkably, while the above j.p.d.f. resembles the j.p.d.f. of 
complex eigenvalues in GinSE (2.3), it is manifestly different from the latter. 

Is it possible to determine the correlation functions for the new complex eigen- 
value model (8.2)? The answer is positive. 

Lemma 8.1. Let Ti. be an n x n random real matrix with no real eigenvalues such 
that its entries are statistically independent random variables picked from a normal 
distribution N(0, 1). Then, the p-point correlation function (1 < p < £) of its complex 
eigenvalues equals 



p("Ho) 

0,p 



(zi,- 



u.Ul r 3 



X pf 



Kt(Zi,Zj) Ki{zi,Zj) 
Ki(Zi,Zj) Ke(z~i,Zj) 



iV2 



exp 



Z2 3+~ Z 1 



2px2p 



Here, n = 2£ and the 'pre-kernel' K£ is 

£-1 

P2j (z)p2j+l (z') - p 2 j (z')p 2 j+l(z) 



£-1 1 

Ke (z,z') = i^2- 



(8.3) 



(8.4) 



The polynomials pj(z) in (8.4) are skew orthogonal in the complex half-plane Imz > 0, 
(P2j+i,P2k)c = -{P2k,P2j+i)c =ir 3 S :j k, (8.5) 

(p2j+l,P2fe+l)c = (p2j,P2k)c = 0, (8.6) 

with respect to the skew product 

(f,9)c = j^^zedc^fj ex p(-^±i!) [f(z)g(z) - f(z)g(z)] . (8.7) 



Proof. By definition (2.7), the p-point correlation function is 

R^\ Zl ,..., Zp ;n) = —^— rj / d 2 z j P 7i0 (z 1 ,---,z e ). (8.8) 

\Z—P)' J=p+ i Jlmz 3 >0 

Since the nx n real matrix with no real eigenvalues has £ pairs of complex conjugated 
eigenvalues, it holds that n — 2£. 

Conceptually, the proof to be presented consists of three parts. First, we concen- 
trate on the j.p.d.f. P-n {zi, ■ ■ ■ , zg) and show that it can be represented in terms of a 
certain quaternion determinant. Second, we prove that the quaternion matrix under 
the quaternion determinant enjoys the projection property (see Definition 2.1). Third, 
we apply the Dyson integration theorem to carry out all (£ — p) integrations in (8.8). 



8 Correlations of complex eigenvalues of a matrix without real eigenvalues 



■58 



Part 1. — As the Vandermonde structure of (8.2) mimics that of GinSE (2.3), it is 
tempting to employ the identity 



n 



\Zi - Zj\ \Zi - z 



Ykzj - Zj ) = det [zf 1 ,^ 1 ] i=li 



,2£ 



(8.9) 



i>j = l j=\ 

that helps us reduce the j.p.d.f. P-h to the form 

2 -n(n-l)/4^ 



(n/2)!n?=ir(j/2) 



det [pi_i(zj)) Pi-ife)] i= i,..., M 



x ^Qerfc 



i\/2 



exp 



7 2 4- ? 2 



(8.10) 



Here, Pfc(z) are arbitrary monic polynomials of degree k. 

The very structure of the matrix under the determinant in (8.10), 



det 



Po(zi) Po(zi) 
Pi(zi) Pi(zi) 



Po(z e ) Po{zi) 
Pi(zi) Pi(zt) 



Pll-l{z\) P2i-2<yZi) Pll-l{zi) Pll-l{zi) 

P2l-\{z\) P2£-l(zi) P2t-l(zt) P2t-l{Z() 

suggests that we introduce a set of quaternions {tp (z), • • • , ip e _ 1 (z)} 7 



^i(z) 



P2j(z)+P2j+l{z) - , P2j{z) -P2j+l(z) 



■ e + 



2i 



ei 



, P2j(z) - P2j + l(z) P2j(z)+p 2 j + l(z) 

+ 2 62 + 2i 63 



(8.11) 



(8.12) 
(8.13) 



whose 2x2 matrix representation reads * 

e[V^)]=[ P2j{ f, P2j{ % 

JWJ [ P2j+l{ z ) P2j+l(z) 

As a result, the above determinant can equivalently be written as 

det \pi-i(zj), Pi-i(zj)} i=1 ... i2e =dct [e[^j_i(^-)]] . =1] 

The latter can be put into a quaternion determinant form (see Corollary 5.1.3 in 
Mehta's book (2004)): 

det [6[</ViM] l=1 ,..., e = qdet [AA] M = qdet [AA] ex£ . (8.14) 

Here, A. is an £ x t quaternion matrix with the entries 

A, = Vi-ite) (8-15) 

* The 2x2 matrices e, are defined as follows: 



1 
1 



ei 



i 
-i 



e2 



1 
-1 



e3 



(? 
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and A. is the dual quaternion matrix whose entries 
are determined by the dual quaternion ^pj(z), 

7 , v P2j(z) +P2 j+1 (z) p 2 j{z)-p2j+l{z) 

W) = 2 60 2l 61 

P2j(z) -p 2 j+l(z) „ P2j(z) +p 2 j+l(z) „ 
■ 62 — e 3 



(8.16) 



such that 



e[^.(^)] = e 2 e[^.(^)] T e 2 - 1 



2i 

i>2j+l(z) ~P2j(z) 
-P2j+i(z) p 2j (z) 



Combining (8.13) and (8.14) into 

det [pi_i(zj), Pi-i(zj)] . =1 _ ^ = qdet [AA 

we obtain 

det tPi_i(«j). Pi-i(«j)]i=i,...,2/ 

Here, the (self-dual) quaternion kernel kt admits a 2 x 2 matrix representation 

-Ke(zi,Zj) -Ke(zi,Zj) 

K t (Zi,Zj) Ke(z l7 Zj) 



(8.17) 
(8.18) 
(8.19) 



* e Yl r j (ldet[k e (z i ,z j )] exr (8.20) 

3=0 



6 [k e (zi,Zj)] = 



(8.21) 



with 



t-i 1 



j=0 ^ 



P2j (z)P2j+l H - P2j («j)P2j+l (z) 



(8.22) 



The set of constants {rj} is not fixed so far. 

The above consideration results in the following expression for the j.p.d.f. P-h 
[Eq. (8.10)] 

2 -n(n-l)/4 n^O r J 

^Wo^ii ■ • • > z e) = — 7Z77^V\ — Tv 1 — 7777^ ^ dct z j)kxt 



x ^Qcrfc 

3=1 



(«/2)! n;=ir(j/2) 



i\/2 



exp 



5.23) 



that serves as a proper starting point for evaluating the p-point correlation function 
R 0p °'(zi, - ■ ■ ,z p ;n) specified in Eq. (8.3). 

Part 2. — Now, we are going to prove that the quaternion kg, satisfies the projec- 
tion property Definition 2.1. To simplify the consideration to follow, we set the so far 
arbitrary monic polynomials Pk{z) to be skew-orthogonal in the complex half-plane, 
Imz > 0, as defined by (8.5) - (8.7). 

Having imposed the skew-orthogonality on Pk(z), with rj € M, we are in the 
position to verify whether or not the projection property for the quaternion ke is 
fulfilled. In accordance with the Definition 2.1, one has to consider the integral 



/ = 



d 2 werk 



'lmw>0 



iV2 



exp 



w 2 + w 2 



@[k e {z 1 ,w)}Q[ke(w,Z2)] (8.24) 
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I = 



that equals 

5e(zi,z 2 ) 6i(zi,z 2 ) 
-Se(zi,z 2 ) -8e(zi,z 2 ) 

Here, the function 8t{z\, z 2 ) is defined by the integral 

' w — w 



(8.25) 



5e(zi,z 2 )= / d z weric 

Jlmw>0 



iy/2 



cxp 



w 2 + w 2 



x [Ki(z 1 ,w)K i (w,z 2 ) - Ke(z 2 ,w)Ke(w,zi)]. (8.26) 

Its evaluation, based on (8.22), (8.7), (8.5) and (8.6), is straightforward, the result 
being 

St{zi, z 2 ) = -Kt{z\, z 2 ) (8.27) 

so that 

I = e[k e ( Zl ,z 2 )]. (8.28) 
Put differently, the quaternion kernel kg satisfies the projection property in the form 



d werfc 



Im w>0 



W — W 

iy/2 



exp 



w 2 + w 2 



Kt{z\,w) Kt(w, z 2 ) = ki{z\,z 2 ). (8.29) 



This is precisely (2.86) of the Dyson integration theorem with a quaternion A = 0. 



Part 3. — The above proof of the projection property for the quaternion kernel ki 
in (8.23) paves the way for carrying out the (t — p) integrations in (8.8). Indeed, 
the integrations therein can be performed by virtue of the Dyson integration theorem 
(Theorem 2.1) since 

6 [k e (zx, z 2 )] d = e 2 9 [ht{zu z 2 )] T e^ 1 = 6 [k i (z 2 , zi)] . (8.30) 
To this end, one has to determine the constant c defined by the integral [see (2.8d)] 

' z — z*~ 



ce 
yielding 



d 2 z erfc 



lmz>0 



iy/2 



exp 



z 2 + z 2 



k e (z,z) 



J 

Jin 



c = I d 2 z erfc 

Im z>0 



iy/2 



cxp 



z 2 + z 2 



Ke(z, z) = I. 



(8.31) 



(8.32) 



The projection property (8.29) combined with the result (8.32) brings the key 
integration identity: 



(f 



1 '■ ,=„4-1 JI > 



j=p+l Jlraz 3>° 



d Zj erfc 



iy/2 



I z 2 + z 2 \ 
exp — 2 —r- J - qdeb[kt(zi,Zj)] txt 



qdet [k((zi, Zj)] pxp 



pf 



Ke(Zi,Zj) Ki{Zi,Zj) 

Ke{zi,Zj) Ke(zi,Zj) 



(8.33) 



Applied to (8.8) and (8.23), it results in 



^0,p (^l> ' ' Z P> "0 — Pn,n T-rn 



u e-i 



3=0 ' 3 



xpf 



UU TU/2) 11 

Kt(Zi,Zj) Ke{Zi,Zj) 
Ki(Zi,Zj) K e (Zi,Zj) 



j^J erfc 



iy/2 



exp 



2px2p 



2? -1- v 2 
6 3 T A } 



(8.34) 



2px2p 



where n = 21. This completes the proof of the Lemma 8.1. 



9 Conclusions and open problems 



61 



8.3. Probability p rh0 to find no real eigenvalues 

The technique exposed in the previous subsection allows us to establish the following 
structural result for the probability to find no real eigenvalues in GinOE spectra. 

Corollary 8.1. The probability p nfi to find no real eigenvalues in spectra of GinOE 
equals 

e-i 

Pnfi = mf^umjl rj - (8 - 35) 

Here, n — 21; the constants rj are defined by (8.5). 



Proof. For n — 21 even, the definition (3.1) translates to 

d 2 z 3 P 7io (z 1 ,---,z l ). (8.36) 



Pn.O = TT / 
„_i Jin 



j=l Jlm Z 3>° 

Given Pn in the form (8.23), the above probability reduces to 

ne-i 
i=o 5 



Pn,n 1 lj=0 r i TT f ,2 , , r- I 

Pnfi = 71 n n p , {I / d Zj qdet [K e (z u Z 

t: llj = l i U/ / J •_ 1 JImz j >0 



x n i erfc ( v £ S") cxp ( v - f[ T £l J- ^ 

Due to the projection property (8.29) of the quaternion kernel Ke, the integration in 
(8.37) can readily be performed [see, e.g., the identity (8.33)] bringing 

n j lmz >o rf2 ^ ^ ^-)i/x/ n crfc (^/t) cxp (- f ^J = ^ ! (8 - 38) 

This establishes the result (8.35) thus completing the proof. ■ 

Remark 8.1. To make the formula (8.35) explicit, one has to know the normali- 
sation constants rj which we have failed to determine explicitly so far. 



9. Conclusions and open problems 

To summarise, an exact formula was derived for the probability p n ,k to find precisely 
k real eigenvalues in the spectrum of an n x n random matrix drawn from GinOE. 
Based on the Pfaffian integration theorem (that can be seen as an extension of the 
Dyson integration theorem to kernels that do not possess the projection property), the 
solution found expresses the probability function p n ^ in terms of zonal polynomials jj. 
This links the integrable structure of GinOE to the theory of symmetric functions 
(Macdonald 1998). Undoubtedly, much more effort is needed to accomplish the 
spectral theory of GinOE. Below, we list some of the open problems that have to 
be addressed. 

(j To the best of our knowledge, this is the first ever random-matrix-thcory observable admitting a 
representation in terms of symmetric functions. 
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(I) Probability Function p nyk and Associated Generating Function G n (z). 
The exact solution for the probability p n _ k expresses the probability in terms of zonal 
polynomials of some complicated arguments. Does a more explicit representation (like 
the one for p n , n -2, see (7.1) and (7.29)) exist for p n ^ A similar question arises in 
the context of the solution for the determinantal generating function G n (z) given by 
(6.11). Can the determinant (6.11) be calculated in a closed form as a function of k, 
n and z? We have a few indications that this is a formidable task. 

(II) Asymptotic Analysis of p, lyk in the Large-™ Limit. The problem we wish 
to pose here concerns the large-n behaviour of p riyk in various scaling limits: (a) 
k <~ n°, (b) n — k <~ n° , and (c) k <~ n 1 / 2 (this scaling is prompted by the result (2.28) 
for an average number of real eigenvalues) . If available, the large-n formulae of this 
kind would facilitate a comparison of our exact theory with existing numerical and 
experimental data, as reported f by Halasz et al (1997) and Kwapien et al (2000). 

(III) Correlation Functions of Complex Eigenvalues andp„. . The solution 
presented for these two spectral characteristics involves specific polynomials which are 
skew orthogonal in the complex plane with respect to a somewhat unusual weight 
function containing the complementary error function. Can those be determined 
explicitly to eventually bring closed formulae for both the correlation functions and 
the probability p„,o? Their large-n analysis would be of great interest, too. 

(IV) Correlation Functions of Both Real and Complex Eigenvalues and 
a Generalised Pfaffian Integration Theorem. The calculation of all partial 
(p, g)-point correlation functions Rp. q k for GinOE matrices with a given number (fc) 
of real eigenvalues (defined as an integral of the j.p.d.f. over all but p real and q 
complex eigenvalues, see (2.7)) is yet another important problem to tackle. Also, 
can the unconditional (p, <7)-point correlation functions R p q = ^2 k R P , q k be explicitly 
determined? We believe that progress in this direction can be achieved through a 
proper extension of the Pfaffian integration theorem: 



Here, the notation of the Theorem 5.1 was used. Notice that one should not assume 
the projection property for the kernel Q n . 

The above list of open problems calls for further research of GinOE that will eventu- 
ally unveil the rich mathematical structures underlying this classical but still largely 
unexplored non-Hermitcan random matrix model. 
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Appendix A. The n x n matrix x 

Al. Symmetries of x & n d useful expansions. — The definition (3.27) of the matrix x 
suggests that it is antisymmetric, 

Xjk = ~Xkj, (A.l) 
and purely imaginary 

Xjk = -Xjk (A.2) 
as confirmed by the formula 

Xjk =i J da{w)e- {w2+ ^' 2 [Re q (w) Im q k (w) + Im qj (w) Re q k (w)] . (A.3) 

Due to (3.16) that relates the skew orthogonal polynomials qj(w) to Hermite 
polynomials Hj(w), the following expansions! are useful for even/odd order Hermite 
polynomials taken at w = x + iy: 

m (— 1V'2 2j ' 

ReH 2m (w) = J2 (2-V (~ 2m ^ y 2j H2m-2j(x), (A.4) 
m_1 c_i y-i2 2 J+ 1 

lmH 2m (w) = J2 ( 2j + i)! (- 2m )^+i y 23+1 H2m-2 3 -i(x), (A.5) 
and 

m , 

, , r„o = v 1: 

(2.7)! 



Rei/ 2m+1 H = £ V , (-2m - l) 2j y 2 > H 2m+1 _ 2j (x), (A.6) 



™ ( , _lV'-i2 2 J+ 1 

Imtf 2m+1 ( W ) = ^ 1 + (-2m - l) 2j+1 ^+ 1 J ff 2m _ 2 ,(x). (A.7) 

Here, (— a)„ is the Pochhammer symbol (a > 0), 

T(a + 1) 



(-a)„ = (-1)" 



r(a-n+ 1)' 



,42. Calculation of the matrix elements X2a,2/3 and X2a+1, 2/3+1- — We claim that 

X2a,2/3 = X2a+1, 2/3+1 = (A. 8) 

for all a = 0, 1, • • • and /? = 0, 1, • • ■. To prove this statement, let us consider x 2 a,2/3- 
In accordance with (A.3), this matrix clement is related to the integral containing 

Reg 2Q (w)Img 2(3 (w) + lmq 2a (w) Re q 2 p(w) ■ 



Having in mind (3.16), (A.4) and (A.5), we observe that the above expression, being 
integrated over the x-part of the integration measure 

da(w) = dx ■ erfc(yV2) %) dy, (A.9) 

nullifies due to the products of Hermite polynomials H even ^ (x) i/odd & (x) of even and 
odd orders. Thus, we conclude that X2a,2/3 — 0. By the same token, the matrix ele- 
ment X2o+i, 2/3+1 is zero as well. 
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A3. Calculation of the matrix elements X2a+i,2/3 and £20,2/3+1 • — We claim that 



X2a+l,2/3 - g 



/>oo 

(2a + 1) L^- 1 (-2 2/ 2 ) + 2 y 2 L^ +1 (-2 y 2 ) 



(A.10) 
(A.11) 



In accordance with the symmetry relation (A.l), it holds 

X2a, 2/3+1 = — X2/3+l,2a- 

To prove the result (A. 10), we start with the definition (A. 3), 

X2«+i,2/3 = i I <faWe- (uj2V)/2 

x [Req 2 a+iM Im 172/3 (w) + lmq 2a+ i(w) Rc q 2 p(w)] , (A. 12) 

where the skew orthogonal polynomials Qj(w) are given by (3.16). Substituting them 
into (A. 12), representing w as w = x + iy and making use of the expansions (A. 4) 
to (A. 7), one can carry out the integration over i£l straightforwardly, due to the 
factorised integration measure (A. 9) and the orthogonality of Hermite polynomials 
on R with respect to the weight exp(— x 2 ). Lengthy but straightforward calculations 
result in 



1 a\ 



X2q+1,2/3 



(la, 13 ~ la-1,0) 



(A.13) 



where J a .f3 is 
7a,/3 = (-If-" 



(2a + 1)^ / dyy 2 ^-^- l e y \v{ C (ys/2) 

(2y 2 V 



2a+l 



■m E - 



j\(2a + l-j)\(j + 2(f3-a)-iy. 



(A.14) 



Notice that for a < —1, the sum in (A.14) is void so that ja<-i.p = 0. Otherwise, 
the series can be summed up in terms of Laguerre polynomials, 

W E .,„_^L_,_, = Ll^-\-2y% (A.15) 



j\(2a + l-jy.(j + 2(P-a)-iy. 
to yield 

l a .^(-lf- a (2a + l)^- / ^^-«)- 1 e«%rfc(^)L^ + T)- 1 (-2 2/ 2 ). (A.16) 

From now on, the Laguerre polynomials of "negative order" are interpreted to be ze- 
ros. Substituting it back to (A.13), we end up with (A. 10). 

Finally, the formula for X2a, 2/3+1 follows from the symmetry relation (A. 11). 

A4- The structure of \. — To summarise the calculations of Appendix A, the structure 
of the n x n matrix x is presented below. For n = 2m even, the matrix x is 



-even 





-X01 



X01 




— Xl,2m-2 

V -X0,2m-1 





Xl,2m-2 



X0,2m-1 




\ 



X2m-2,2m-l 
_ X2m-2,2m-l 



(A.17) 
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For n = 2m + 1 odd, the structure of x follows the pattern 



y° dd 





-Xoi 



"X0,2m-1 




Xoi 






— Xl,2m 



Xo,2m-l 





o \ 

Xl,2m 



X2m-l,2m 
-X2m-l,2m j 



(A.18) 



Appendix B. The n x n matrix a = 2i (fix) 

Bl. Calculation of the entries cr a ,f3 = 2i (/tx)a,/3 for n even. — The matrices ft, x and 
thus (fix) are sensitive to the parity of n, the matrix dimensionality. For n — 2m, we 
make use of (3.21), (A.8), (A.ll) and (A.17) to derive 

(Ax)& + i = (Ax) e 2^i, 2/3 = o, (B.l) 

(Ax)l™2/3 = -T- X2a+l,2/3 = T X2/3,2a + l, (B.2) 
(Ax)2a+l,2/3+l = T- X2a, 2/3+1 = T X2/3+l,2a = (Ax)2™2a- ( B - 3 ) 

Here, a = 0, 1, • • • and /3 = 0, 1, • • •. 

An integral representation for the entries <72™2/3 an d ^20+1,2/3+1 can be read off 
from (A. 10). Explicitly, we have 

poo 

~lf- a / dyyW-^- 1 e y2 er{c(yV2) 
Jo 

x "(2a + 1) Ll^-\-2y*) + 2 y 2 L^ +1 (-2 y 2 )] (B.4) 

-even ~cven /T3 r\ 

°2a+l, 2/3+1 — T~ "213. 2a- \ a -°) 

All other entries are zeros. 

Yet another representation in terms of 7 aj/ g given by (A. 16) is useful. Equations 
(A.13), (B.2) and (B.3) yield: 

^2/3 = ^ (7a,/3 ~ 7a-l,/j) , (B.6) 

(3\ 

^2a+l,2/3+l = J" (7/3,a - 7/3-l,a) • (B.7) 

All other entries nullify in accordance with (B.l). 

B2. Calculation of the entries a a jj — 2i (fix)a,f3 for n odd. — For n = 2m + 1, we 
make use of (3.22), (A.8), (A.ll), (A.18), (B.2) and (B.3) to derive 

(Ax)2tVi = (Ax)r+i, 2/ 3 = 0, (B.8) 

I m— 1 , 

= (1 - ^, m )(Ax)S - F 1 S ">™ E if (AX)^' ( B -9) 



; 2a, 2/3 



and 



(Ax) 



odd _ h P 

2a+l, 2/3+1 - 7— 



777,! 

(Ax)2™2a — 7; (AX)I^277 



(B.10) 
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Here, a. = 0, 1, • • • and (3 = 0, 1, • • •. 

The above formulae simplify if expressed in terms of 7 aj/ g akin to (B.6) and (B.7). 
Straightforward calculations lead to the following result for the entries of <r odd : 

a! r . i 

(1 - £>a,m)(7a,/3 _ 7a-!,/?) _ &a,m lm-1,/3 , 



-odd 



'2a,2/3 - T- 
-odd 

7 2a+l, 2/3+1 



T-(7/3,a - 7/3-1, a) - J— (7/3,m ~ 7/3-l,m) 



(B.11) 
(B.12) 



Finally, the integral representations for (B.ll) and (B.12) can be obtained from (A. 16). 

Appendix C. The [n/2\ x [n/2\ matrix g and the trace identity 

CI. The definition of g. — Since the n x n matrix <x has half of its entries vanishing, 
it is useful to define a reduced, [n/2\ x [n/2\ matrix g, such that 



-even / -i \8— a -even / i \/3— a /~ ~ \ 

Q a ,(i = (-1) CT 2«,2,3 = (-1) — (7a,/3 - 7a-l,/3) 



for n = 2m even, and 



-odd / i\f3-ahf3 -odd / i\fl-Q 

ffa,/? = T ^2/3+1.20+1 = M) 

"a 

_ -even i \/3-m "Zl ^ 



T 2a,2/3 



m! /i^ 



(CI) 



even 



h m [31 



'Jr. 



(C.2) 



for n = 2m + 1 odd. Explicitly, 



lyW-^-^v 2 erfc(y\/2) 

-2(/3-a)-l 



(2a + 1) Lr + "i a;_1 (-2l/ a ) + 2y 2 ^<?_r )+1 (- V) (C.3) 



and 



-odd -even 

Qa,f3 — Qa,f3 



k ^ (2m)! /? 



| Qa.,m- 



(C.4) 



The above formulae have been obtained from the definitions (C.l) and (C.2) with the 
help of (B.4) and (3.18). 

C2. The trace identity. — It turns out that the probability function p n ^ is most eco- 
nomically expressed in terms of g. To realise this, it is instructive to prove the following 
trace identity. 

Lemma CI. Let & be an n x n matrix with the entries given by either (B.6), (B.7) 
or (B.ll), (B.12) depending on the parity of n. Then, the trace identity 

tr (o.n-i)^ = 2tr (0iLrl/ /2j_i)£? : ' (C.5) 
holds for all j = 1, 2, • • • and the matrix g defined by (C.l) and (C.2). 

Proof. Since the matrix elements <7 2aj 2j a+i and <t 2q ,+i.2/3 are zeros, the trace of a 3 
(j = 1, 2, • • •) can always be separated into two pieces: 



tr (0,n-l)<T j - tr (0 J„/ 2 ]-l)a J +tr( ,[n/2j-l) b • 



(C6) 
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The matrices a and b are denned through their entries, 

a a ,f) = &2a, 2/3, (C.7) 
bajj = &2a+l, 2/3+1, (C-8) 

and explicitly depend on the parity of n in accordance with the discussion in Appendix 
B. 



(i) The case n — 2m even is the simplest one because of the relation (B.5). Indeed, 
its straightforward use results in 



tr 



(0,m- 



7 'J 



This can further be simplified due to(C9), (C.7), (C.6) and (C.l), 



tr (0,2m-l)O- J = 2 t r (o,m-i)a J 



(C.9) 
(CIO) 



This completes the proof of (C.5) for n = 2m even. 



(ii) The case n = 2m + 1 odd is a bit more complicated since a simple analogue 
of (B.5) does not exist. Instead, we have [(C.7) and (C.8)] 

a a ,p — C a (1 — C>a,m)(7a,/3 — 7a-l,/3) — #a,m7ro-l,/3 , (Cll) 
b a ,0 = C a (73^ - 7/3-1, a) - Cm {l0,m - 7/3-l,m), (C12) 

where c Q = a\/h a has been defined in (3.23). In writing (C.12), we dropped the 
prefactor (5\/a\ appearing in (B.12) since it does not affect the value of the second 
trace in (C.6). To prove (C.5), we will start with (C.6) in order to demonstrate that 



tr (0,m-l)b J 



(C.13) 

This will be followed by a proof that either of these traces reduces to tr (o,m-i)£° dd J - 
Let us prove (C.13) by focusing on the eigenvalues {A} of the matrices a and b 
which are the roots of the secular equations 

det[a-Al] (m+1)x(m+1) =0, (C.14) 
det \b- Ail =0. (C.15) 

L Jmxm v ' 

We claim that (C.13) holds because exactly one out of (m + 1) eigenvalues of the 
matrix a is zero whilst the remaining m eigenvalues of a coincide with m eigenvalues 
of b. Put differently, we are going to prove that 

det [a -Ail, ^ , . ^ = — A det [6 — A 11 . (C.16) 

L J (ro+l)x(ro+l) L Jmxm v 1 

The proof consists of four steps. 

Step \ 1. Consider the (m + 1) x (m + 1) matrix a — Ai under the determinant 
in the secular equation (C.14), 



/ Co 7o,o - A 

ci ri i0 



co7o,i 
ci Ti t i — A 



i r rn — 1 5 o c m _i r rn 

Cm / 77T, — 1,0 ~~ Cm im —1,1 



70,m-l 

ci I\ m _i 

-l r m _i_ m _i a 

"Cm ] m— l,m— 1 



CO 70,n 

ci i\ r 



\ 



Cm— i r rn _ i ; m 

~ Cm 7m — 1 , m / 



(C.17) 
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where T at p = 7 Q)J g — 7 a _i j; g. Let us perform a number of operations with rows and 
columns that will leave the value of the secular determinant intact. First, we multiply 
the content of the first row by Ci /co and add it to the second row. Having done this, 
we multiply a new content of the second row by C2/C1 and add it to the content of 
the third row. We go on with this procedure until we arrive at the modified mth row 
whose content is multiplied by c m /c m _i and further added to the last, (m + l)th, 
row. While not affecting a value of the determinant in (C.14), the above sequence of 
transformations brings (C.17) to the form 

/ co 7o,o -A c 70,1 ••• c 7o,m \ 

ci 71,0 - (ci/c )A ci7i,i-A ••• ci7i, m 



Cm-l 7m-l,0 - (c m -l/co)A Cm-iTVn-i,! - (c TO _i/ci)A 

-(c m /c )A -(c m /ci)A 



Cm — 1 7m-l,m 

-(c TO /c m )A J 

(C.18) 

The last row of this equation suggests that a factor A can be taken out of the secular 
determinant. In other words, A = is always an eigenvalue of the (m + 1) x (m + 1) 
matrix a. 



Step t| 2. Next, we multiply the content of the first row in (C.18) by ci/c and subtract 
it from the content of the second row; having done that, we multiply the content of 
the modified second row by C2/C1 and subtract it from the third row; going on with 
this procedure, we arrive at the (m — l)th row, multiply it by c m _i/c m _ 2 to subtract 
this from the content of the mth row. We do not touch the last, (m + l)th row. This 
set of transformations yields 

/ c 70,0 -A c 70,1 ••• C 7o,m-l C 7o,m \ 

ci Ti.o ci Ti 1 — A ••• ciri. m _i 



Cm-l Tm-1,0 C m _i T m _ii 

V -(c m /c )A -(c TO /ci)A 



CO 70,m 
Cl Ti m 



— (c m /c m _i)A — (c m /c m )A J 

(C.19) 

Note, that all but the last row of the matrix (C.19) coincide with those in (C.17). 



Step t]3. Now, let us factor out c from the first row, c\ from the second row, 
c m _i from the mth row to obtain 

m— 1 

n ^ 

3=0 



times 

/ 7o,o - A/c 70,1 

I\o r M -A/ci 

lm-1,0 Tm-1,1 

V -(c TO /c )A -(c m /ci)A 



70,m-l 
lYm-l 

Tm— l,m— 1 A/c m _i 

— (c m / c m -i)A 



70,m 
Tl,™ 

^m- l,m 

-(c m /c m )A J 



(C.20) 
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Next, we multiply the first column by Co, the second column by c\, the mth column 
by c m _i, and do not alter the last, (m + l)th column. This leads us to 

/ c 7o,o-A Ci7o,i 



co r m -i,o 

\ — c m A 



— c m A 



c m-l 70,m-l 
Cm— 1 Tl-m — 1 

— 1 Tm— 1 

— c m A 



70,m \ 

ri. m 



-A 



(C.21) 



Step t] 4. Now, we multiply the last, (m + l)th column by c m and subtract it from the 
first, second, the mth column to reduce (C.21) to 



/ C 7o : - Cml0,m ~ A 
CO IYo — C m Ti m 



Co r m — i,o c m r m _i_ m 
o 



70, m \ 
Tl,m 



1 m — l.m 

-A / 



(C.22) 



The determinant of the latter matrix can be calculated via expanding with respect to 
its last row. As a result, the secular equation (C.14) is reduced to 

-A • dct [c/3 (ja.p - Ja-1,0) - C m (j a ^ m - %-l, m ) ~ A S a ,p] ( Q ,/3)=o,- . .,m- 1 • 

A comparison with (C.12) allows us to rewrite (C.23) in the form 



-A dct 

This establishes (C.13). 



Al 



(C.23) 
(C.24) 



Finally, it remains to show that, for all j — 1, 2, • • •, the identity 



tr 



(0,771-1) 



b J =tr 



(0,1 



X odd j 

-i)Q 



(C.25) 



holds. That (C.25) is indeed true, follows from (C.8) and (C.2). This completes our 
proof of the Lemma. ■ 



Appendix D. Calculation of the trace tr(o ^n/2j-i)£ 

Since the matrix g is sensitive to the parity of n, two separate calculations are needed. 

(i) The case n = 2m even. — To calculate the trace, we make use of (C.3) to write 
down 

m—l 



tr 



(0,771-1)5 



^ ] Qa.a 



dyy 1 e v evfc(yV2) 



o 



Q=0 

m—l 
a=0 

First, to put (D.l) into a more tractable form, we apply the identity § 



Ln m H = 



w" 



{-n) r 



(D.l) 



(D.2) 
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which, in the context of (D.l), reads 

^ lW = "&TT^ W ' (D ' 3) 

The use of (D.3) reduces the integrand of (D.l) to 

(2a + 1) L 2Q 1 +1 (-2y 2 ) + 2y 2 l} 2a _ x (-2y 2 ) = 2y 2 [L\ a {-2y 2 ) + L^^y 2 )] . (D.4) 
Second, we spot that the transformation 

L x v -\w) = Ll{w)-Ll_ l {w) (D.5) 
applied to (D.4), yields 

L\ a (-2y 2 ) + LL_ 1 (-2y 2 ) = L 2 2a (-2y 2 ) L 2 2a _ 2 (-2y 2 ). (D.6) 



T 2 7-2 

L „-L "2a-l _ ~2a-2 



As a consequence, the summation in (D.l) can be performed explicitly, 

m— 1 

]T [(2a + 1) i 2Q 1 +1 (-2y 2 ) + 2y 2 L\ a _ x {-2y 2 ) 

m—1 

= 2y 2 J2[LL(-2y 2 )-L 2 2a _ 2 (-2y 2 )_ 



a=0 



a=0 

= 2y 2 L 2 2m _ 2 (-2y 2 ), (D.7) 
resulting in a remarkably simple formula 

tr (0 , m _i)e° vcn -2 / dyyey 2 eric(yV2)L 2 2m _ 2 (-2y 2 ). (D.8) 
Jo 



(ii) The case n = 2m + 1 odd. — To calculate the trace, we combine (C.l) and (C.4) 
into 

Qa,a = Qa.a ~ T {la,m ~ 7a-l,m) ■ (D-9) 

Summing it up, we derive 

4_ ~odd i ~even - /t-\ 1 n\ 

t r (O,m-l)0 = tr (O,m-l)0 - T 7m-l,m- (D.10) 

Further use of (D.8) and (A. 16) yields 

/•GO 

tr (0 . m -i)e° dd = 2 dyyey 2 erfcQ/v^) [iL- 2 (-2y 2 ) + L^^y 2 )] . (D.ll) 
Jo 

With the help of the identity (D.5), this eventually simplifies to 

p OO 

tr ( o, m -i)^ odd = 2 / dyye ! ' 2 crfc(yv / 2)LL-i(-2y 2 ). (D.12) 
Jo 



The formulae (D.8) and (D.12) can be unified into a single equation, holding for 
n of arbitrary parity: 

P OO 

tr (0 , L n/2j-i)e = 2 / dyy e ^crfc(j/\/2)i 2 _ 2 (-2y 2 ). (D.13) 
Jo 
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